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This  thesis  develops  a  straightforward  and  accurate 
method  to  determine  mechanical  heating  system  changeover 
timing.  Transitional  season  weather  causes  buildings  to 
stay  within  acceptable  thermal  comfort  ranges  without 
requiring  supplemental  heating.  Significant  energy  savings 
can  be  realized  if  heating  plants  can  be  shut  down  during 
periods  they  are  not  required.  The  term  "trigger 
temperature"  is  used  to  describe  the  minimum  acceptable 
comfort  conditions  which  must  be  maintained  before 

I  • 

changeover  can  be  made.  Computer  simulation  of  ninety-four 
hypothetical  office  buildings  is  used  to  record  data  for 
hourly  building  temperatures.  Statistical  regression 
analysis  techniques  are  employed  to  analyze  the  data  and 
develop  temperature  prediction  models  for  a  range  of 
building  envelopes  and  occupancy  schedules.  Hand 

calculation  and  micro-computer  applications  are  offered  to 
provide  changeover  decision  criteria  for  facility  managers 
using  this  method. 

Key  words  fo^'this  thesis  are:  transitional  seasons, 

changeover  timing,  regression  analysis,  thermal  comfort,  and 

« 
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CHAPTER  ONE 


INTRODUCTION 


Since  Che  Oil  Embargo  of  1973,  social,  economic,  and 
pollcical  pressures  have  underscored  the  need  for  energy 
conservation  throughout  the  world.  The  recent  reduction  in 
crude  oil  prices  has  been  directly  linked  with  the 
conservation  efforts  at  work  during  the  past  year.  The 
resultant  reduction  In  energy  demand  has  created  an  oil 
surplus:  a  gap  between  the  demand  and  available  supplies. 
One  major  contribution  in  reducing  energy  consumption  has 
come  from  the  energy  conservation  accomplishments  In 
building  design  and  operation.  From  the  Federal  Energy 
Administration  Indoor  temperature  guidelines  established  In 
1974,  there  has  been  a  conscientious  effort  to  make 
buildings  more  "energy  efficient"  without  sacrificing  the 
comfort  of  the  occupants.  Increasing  attention  Is  being 
given  to  annual  energy  costs  of  existing  and  proposed 
facilities . 

Many  different  efforts  have  been  directed  toward 
reducing  building  energy  consumption.  Improvements  have 
been  made  in  the  operating  efficiency  of  heating, 
ventilation,  and  alr-cond itlonlng  (HVAC)  equipment,  and 
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facility  managers  and  owners  have  concentrated  on  tightening 
the  building  thermal  envelope.  Engineering  analysis  methods 
have  been  developed  to  maximize  energy  efficiency.  These 
techniques  range  from  simple  energy  audits  which  track  fuel 
consumption  to  sophisticated  computer  controlled  energy 
monitoring  systems.  One  area  which  has  not  received  much 
attention  is  that  related  to  facility  energy  use  during 
transitional  seasons. 


Background 

A  critical  weather  transition  period  occurs  each  Spring 
and  Fall  for  facilities  in  temperate  climates.  At  this 
time,  the  startup  or  shutdown  of  the  facility  heating  plant 
must  be  considered.  The  timing  of  this  "changeover”  is 
critical  because  it  directly  affects  the  annual  energy 
consumption  of  the  building.  In  addition,  it  is  often  an 
irreversible  decision  which  may  adversely  limit  the  comfort 
for  the  facility  occupants.  It  is  Important  that  the  timing 
of  the  changeover  reflects  a  careful  balance  between 
maintaining  an  acceptable  thermal  comfort  within  the 
building  and  the  economic  pressures  to  save  on  energy  bills 
when  fuel  costs  are  high.  The  potential  for  savings  lies  in 
a  changeover  timing  which  avoids  costs  associated  with 
running  and  maintaining  a  heating  plant  during  a  period 
without  any  heating  requirements.  When  the  facility  is 
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heated  by  a  central  plant,  the  savings  potential  Is  the 
greatest . 

Where  a  central  heating  plant  serves  several 
facilities,  the  changeover  timing  has  additional  Impact  on 
energy  use  and  potential  savings.  Central  plants  must  be 
prepared  for  heating  load  requirements  which  may  not  occur. 
In  addition,  losses  occur  In  the  distribution  network 
serving  all  facilities.  Central  heating  plants  are  often 
constructed  where  one  organization  owns  and  operates  all  of 
the  facilities  In  the  distribution  network.  This  situation 
Is  common  on  university  campuses  and  military  Installations 
as  well  as  In  Industrial  complexes,  yet  no  simple  method  has 
been  created  for  changeover  estimation.  For  Instance,  the 
United  States  Air  Force  has  no  existing  method  to  use  for 
Its  numerous  (over  100)  locations  around  the  world.  Nor 
does  ASHRAE,  which  provides  many  of  the  engineering  design 
standards  for  HVAC  equipment,  describe  a  changeover 
estimation  method  In  any  of  Its  publications.  The  purpose 
of  this  thesis  Is  to  propose  a  procedure  from  which  a 
changeover  decision  process  may  be  developed. 

Objectives 

The  objective  of  this  thesis  is  to  develop  a  simple  and 
accurate  method  to  determine  changeover  dates  for  office 
type  buildings,  by  predicting  how  the  buildings  will  perform 
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during  transitional  seasons.  It  Is  Intended  that  this 
method  would  be  used  by  facility  manager  and  engineer  alike 
when  making  building  energy  analyses. 

Research  will  be  comprised  of  the  following  steps: 

1.  Evaluate  facility  manager  needs  for  a  changeover 
methodology  from  Input  information  to  end  use. 

2.  Review  existing  changeover  estimation  techniques 
and  their  present  applications. 

3.  Use  heat  balance  engineering  analysis  methods  to 
determine  how  a  range  of  office  type  buildings  perform 
during  the  transitional  seasons. 

4.  Determine  the  moat  critical  weather,  building  en¬ 
velope,  and  occupant  use  characteristics  affecting  the 
building  changeover  conditions. 

5.  Employ  statistical  regression  techniques  to  predict 
Internal  building  temperatures. 

6.  Integrate  weather  and  user  Input  data  into  the 
model  developed  and  establish  Its  accuracy  and  range  of 
application. 


Limitations 

The  methodology  developed  in  this  thesis  is  limited  to 
determining  the  changeover  date  for  mechanical  heating 
systems  only.  The  changeover  decision  Is  based  upon 


building  temperatures  during  the  hours  of  6am  to  10am. 
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Occupancy  schedules  used  In  this  research  have  been 
developed  by  studying  the  weekdays  of  Tuesday  through  Friday 
In  a  five-day  work  week.  Data  generated  for  the  statistical 
portion  of  the  research  are  limited  to  slab  on  grade  single 
story  medium  density  masonry  structures  with  office  type 
occupancies.  Outdoor  temperatures  used  have  dally  ranges  of 
from  10  F  to  30  F,  and  solar  data  Is  typical  of  from  35  to 
43  degrees  North  or  South  latitude  locations.  Applications 
for  weather  or  building  types  outside  these  limits  have  not 
been  validated  by  the  statistical  model  developed  as  part  of 
this  thesis. 


Contribution 

The  contribution  of  this  thesis  Is  the  development  of  a 
method  to  quickly  and  accurately  determine  building  heating 
system  changeover  dates.  This  methodology  also  provides  the 
user  with  a  technique  to  assess  the  impact  of  differing 
changeover  dates  on  each  facility  by  predicting  how  each 
will  perform  without  a  heating  plant.  Although  established 
and  accepted  methods  are  available  to  calculate  building 
performance  under  these  conditions,  there  Is  no  simple  yet 
accurate  approach  which  accounts  for  heat  gains  due  to  solar 
radiation  and  the  time  lag  effects  of  building  thermal  mass. 
These  thermal  gains  and  losses  are  critical  If  an  accurate 
building  temperature  Is  to  be  calculated  for  use  In  making 
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changeover  decisions.  The  development  of  a  simple  and 
precise  model  to  determine  building  changeover  dates  offers 
the  area  of  energy  management  a  valuable  tool  to  conserve 
energy  without  sacrificing  building  occupant  comfort. 

Overview 

The  development  of  this  thesis  Is  presented  In  the 
following  chapters:  Chapter  Two  reviews  facility  manager 
needs  when  estimating  building  changeover  and  describes 
possible  engineering  applications.  Chapter  Three  reviews 
several  existing  methods  which  could  be  used  to  make 
changeover  date  estimations.  Chapter  Four  details  the 
statistical  analysis  procedures  which  are  used  to  establish 
building  temperature  prediction  models.  Chapter  Five 
Integrates  the  statistical  models  Into  the  Trigger 
Temperature  Method  for  determining  changeover  timing  and 
provides  examples  to  demonstrate  the  method's  capabilities 
and  applications.  Chapter  Six  offers  conclusions  and 
recommendations  for  future  research. 
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CHAPTER  TWO 

FACILITY  MANAGER  NEEDS 

The  responsibility  for  an  organization's  building 
operation  and  maintenance  program  Is  typically  under  the 
central  control  of  the  facility  manager.  The  facility 
manager  may  direct  the  needs  for  a  single  building  or  an 
entire  Industrial  complex.  In  order  to  maximize  potential 
energy  savings,  a  changeover  method  should  provide  him  with 
a  means  to  ensure  buildings  are  heated  adequately,  as  well 
as  efficiently.  Energy  savings  can  be  realized  If  heating 
plants  can  be  shut  down  during  periods  they  are  not 
required.  An  Ideal  analysis  method  would  offer  flexibility 
in  approach,  be  straightforward  and  simple  In  application, 
yet  have  a  high  degree  of  accuracy  when  estimating 
changeover  timing. 

Method  Flexibility 

Changeover  concepts  are  addressed  by  large  and  small 
organizations  alike.  Although  the  scope  of  application  may 
differ,  any  changeover  method  would  need  to  Include 
buildings  of  varying  size,  construction,  occupant  use,  and 
site  conditions.  It  should  apply  to  older,  less  efficient 
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buildings,  as  well  as  newer,  energy  efficient  ones.  As 
older  buildings  are  upgraded  to  be  thermally  more  efficient, 
their  performance  during  the  changeover  seasons  Is  altered. 
A  system  of  changeover  analysis  should  be  able  to  adjust  to 
these  energy  upgrades. 

Application  Simplicity 

A  changeover  method  should  be  straightforward  and 
easily  applied  for  a  multitude  of  facilities.  Facility 
management  organizations  are  not  always  staffed  to  do 
detailed  energy  analyses  and  must  concentrate  primarily  on 
the  maintenance  and  running  of  the  "physical  plant”.  This 
puts  a  limit  on  the  time  and  engineering  expertise  which  may 
be  available  for  changeover  estimation.  Each  hour  added  to 
this  estimation  process  by  extensive  and  time  consuming 
engineering  calculations  Is  an  hour  which  can  be  spent  In 
other  areas  of  energy  conservation. 

Typically,  facility  management  for  smaller 
organizations  Is  not  staffed  with  engineers  to  perform 
changeover  analysis.  This  responsibility  may  rest  with 
central  heating  plant  managerc  or  individual  building 
custodians.  They  need  a  method  which  does  not  rely  upon  the 
expertise  of  architects  or  engineers  trained  In  energy 
analysis  techniques.  Simplicity  and  ease  of  use  Is  their 
primary  concern.  Regardless  of  the  end  user  within  a 
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facility  management  organization,  a  changeover  method  should 
offer  an  accuracy  level  sufficient  to  ensure  thermal  comfort 
for  the  building  occupants. 

Accuracy  and  Accountability 

When  estimating  changeover  timing,  the  accuracy  of  any 
method  la  a  function  of  several  parameters.  The  most 
accurate  method  of  analysis  would  be  one  derived  from 
observations  of  actual  facilities  In  operation.  This  is  an 
Impracticable  approach,  however,  because  it  would  not  be 
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used.  For  example,  weather  data  Is  one  uncontrollable 
factor  which  can  limit  the  accuracy  of  an  otherwise  sound 
analysis  method.  Another  accuracy  constraint  is  the  time 
Interval  for  which  a  method  predicts  changeover. 

A  method  which  determines  changeover  dates  on  a  monthly 
basis  does  not  offer  the  sensitivity  needed  to  be  of  much 
value  to  the  user.  On  the  other  hand,  one  which  can 
pinpoint  this  date  to  a  single  day  of  the  month  may  not  be 
considering  dally  weather  fluctuations.  The  Ideal  method 
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would  consider  at  least  weekly  weather  trends,  yet  retain 
sufficient  precision  to  determine  occupant  comfort 
conditions . 

Occupant  comfort  is  the  most  important  accountable  item 
in  a  changeover  analysis.  The  concept  of  a  changeover  date 
negates  all  of  the  reasons  buildings  are  heated  in  the  first 
place  if  it  cannot  in  some  way  estimate  whether  or  not 
people  will  be  comfortable.  Occupant  comfort  should  be 
recorded  in  some  way  so  that  changeover  conditions  can  be 
determined  as  a  function  of  a  specifed  comfort  condition. 
There  are  several  methods  which  can  be  used  to  assess 
thermal  comfort  for  occupants.  The  relative  merits  for 
these  methods  are  covered  in  Chapter  Three. 
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CHAPTER  THREE 

CHANGEOVER  ESTIMATION  TECHNIQUES 

Building  changeover  dates  can  be  estimated  In  several 
ways.  The  historical  data  approach  uses  records  of 

changeovers  from  previous  years  as  a  data  base  from  which 
future  predictions  may  be  made.  Engineering  calculations 
may  be  made  for  either  steady  state  or  transient  thermal 
conditions  to  determine  the  optimal  changeover  date.  The 
end  goal  in  all  cases  is  to  assess  when  weather  conditions 
have  established  a  situation  where  there  Is  no  longer  a 
facility  heating  requirement. 

Historical  Data 

The  simplest  approach  uses  historical  data  as  the  basis 
for  estimating  changeover  timing.  This  general  procedure 
employs  past  records  of  weather  and  building  performance  to 
predict  future  trends.  The  number  of  years  observations 
have  been  taken  is  an  important  facet  of  this  approach. 
Records  for  each  facility  are  kept  of  past  years'  changeover 
dates.  When  the  heating  plant  Is  turned  off  (or  on),  a 
record  Is  made  of  how  this  timing  affected  the  facility 
occupants.  The  facility  manager  assumes  that  if  no 
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complaints  are  received,  the  changeover  date  chosen  was 
acceptable.  Weather  data  for  this  period  is  recorded  and 
compared  with  weather  from  previous  observations.  The  goal 
for  the  coming  seasonal  changeover  Is  then  selected  as  a 
compromise  between  energy  conservation  goals  and  the  number 
of  "too  cold”  complaints  from  the  occupants. 

In  application,  this  management  by  exception  approach 
has  several  limitations.  Building  temperatures  cannot  be 
determined,  and  occupant  thermal  discomfort  can  be  predicted 
only  Indirectly.  In  addition,  accuracy  in  this  approach  is 
primarily  a  function  of  the  quantity  of  previously  recorded 
observations.  Unfortunately,  this  method  rarely  accumulates 
enough  sample  years  to  be  of  much  value. 

As  a  facility's  use  changes  or  architectural 
modifications  are  made,  the  thermal  performance  of  the 
facility  also  changes.  Added  insulation  or  new  occupancy 
hours  will  change  how  a  building  will  react  to  the 
environment.  As  a  result  of  any  changes,  the  previous 
years'  records  become  invalid  predictors  for  future 
performance.  With  recent  efforts  to  insulate  and  improve 
thermal  performance  of  his  buildings,  a  facility  manager  may 
find  that  it  is  never  possible  to  accumulate  enough  data 
years  using  the  historical  data  method.  A  more  useful  and 
scientific  approach  Involves  using  engineering  heating  load 
calculations  to  estimate  changeover  timing. 


Steady  State  Analysis 


Steady  state  analysis  is  the  most  widely  used  and 
accepted  method  to  determine  heating  transfer  rates  (ASHRAE, 
1981).  Steady  state  building  load  calculations  are  made  by 
assuming  that  the  rate  of  heat  flow  through  a  building 
component  Is  constant  over  a  period  of  time  and  that  the 
rate  of  flow  can  be  approximated  by  comparing  the  potential 
energy  differences  at  the  beginning  and  end  of  the  period 
(usually  one  hour).  The  thermal  transmittance  ,  or  U-value, 
is  the  coefficient  of  heat  transmission  for  the  material 
through  which  the  heat  Is  passing  (ASHRAE,  1981).  Table  3-1 
lists  ASHRAE  steady  state  equations  used  for  estimating 
building  heat  gains  and  losses. 

Balance  Point  Temperature 

The  weather  condition  which  causes  a  building's  heat 
gains  to  be  equal  to  its  heat  losses  is  referred  to  as  the 
facility's  balance  point.  The  outside  dry  bulb  temperature 
under  these  conditions  is  called  the  balance  point 
temperature  (BPT).  For  residential  construction,  this  has 
typically  been  65  F,  but  the  greater  internal  loads  due  to 
lighting  and  occupants  cause  office  type  facilities  to  have 
much  lower  values  (ASHRAE,  1981).  By  calculating  the 
balance  point  temperature,  conditions  can  be  determined  when 
no  net  heat  will  be  required. 
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Table  3-1 

Steady  state  heat  transfer  equations 


Component 
Envelope 
Infiltration 
Ventilation 
Occupanc  y 

Lights  and  equipment 


Steady  state  equation 
U*A*( IDBT-ODBT) 

1 . 1*CFM*( IDBT-ODBT) 

1 . 1*CFM*( IDBT-ODBT) 
255*PEOPLE 
3.413*WATTS 


Units 

Btuh  (  -  ) 
Btuh  (  -  ) 
Btuh  (  -  ) 
Btuh  (  +  ) 
Btuh  (  +  ) 


where  *  denotes  multiplication 

Btuh  =  Btu/hour  heat  gains  (  +  )  or  losses  (  -  ) 
UA  «  U-value*area  of  the  envelope  component 
IDBT  •  Indoor  dry  bulb  temperature 
ODBT  «  outdoor  dry  bulb  temperature 
PEOPLE  ■  number  of  facility  occupants 
WATTS  ■  the  total  Internal  electric  load  In  watts 


Note  ;  255  Btuh  per  person  is  sensible  heat  only. 
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An  application  of  the  BPT  concept  can  be  demonstrated  In  an 
example  problem: 

At  the  balance  point  temperature,  building  heat  losses 
are  matched  by  gains,  where: 

Heat  losses  «  Transmission  +  Infiltration 
Heat  gains  ■  Lights  +  People  +  Solar 

Table  3-2  lists  the  building  data  for  a  2000  sf  modern 
construction  office  building. 


Table  3-2 

Balance  Point  example  building 


Component 


Value 


Transmission  Area, 

Walls  1600 

Glass  200 

Roof  2000 


sf  U-Value,  Btu/ ( hr* s f * F) 
0.10 
0.55 
0.045 


Infiltration 

Occupants 

Lights 


130  cu  ft/mln 
16  people 
2.2  watts/sf 


Using  this  data  to  calculate  the  building's  balance  point 
temperature,  the  building  Indoor  temperature  must  first  be 
set.  If  an  Indoor  setting  of  72  F  Is  used. 

Heat  losses  In  MBtuh  ■  50 3*( 7 2-ODBT ) 

Heat  gains  In  MBtuh  19,097.0 
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By  setting  the  gains  and  losses  equal, 

ODBT  »  balance  point  temperature  •  34.0  F  -  Solar 

The  balance  point  method  has  some  limitations  when  used 
for  estimating  changeover  timing.  The  thermal  gains  from 
solar  and  the  time  delayed  gains  from  the  building  mass 
cannot  be  determined  by  steady  state  analysis  techniques. 

Solar  Gains 

Thermal  gains  due  to  the  sun  can  be  a  significant 
portion  of  the  total  heat  gain  for  a  facility.  Solar  gain 
alone  may  reduce  the  annual  heating  load  by  10  -  15  !?  for  an 
office  building  having  an  elght-to-one  ratio  of  building 
floor  area  to  south^faclng  glass  area  (Balcomb,  1980). 
However,  due  to  the  transient  nature  of  solar  conditions, 
these  gains  cannot  be  calculated  using  stead/  state 
equations.  Figure  3-1  shows  how  the  load  profiles  for  the 
previous  example  problem  would  be  affected  If  solar  gains 
were  considered.  Steady  state  analysis  methods  are  also 
unable  to  account  for  the  thermal  capacitance  effects  caused 
by  a  building's  mass. 

Thermal  Capacitance  Effects 

The  capacitive  effect  of  a  building  Is  another 
significant  contribution  to  the  building  thermal  balance. 
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LIGHTS  OH  - ;  LIGHTS  OFF 


Figure  3-2.  Thermal  mass  storage  effect 


Each  hour  of  the  day,  thermal  energy  is  being  stored  and 
released  by  the  building  roof,  slab,  walls  and  partitions. 
During  the  occupied  hours,  this  building  mass  absorbs  heat 
from  the  lights  and  sun  as  well  as  the  ambient  air.  After 
the  lights  are  turned  off  the  stored  heat  is  released  to  the 
space  as  the  processes  reverses  during  the  unoccupied  hours. 
The  amount  of  heat  which  can  be  stored  In  and  released  by 
the  building  thermal  mass  will  depend  upon  many  factors, 
Including  the  building  construction  and  Interior  layout. 

Though  widely  used  for  load  calculations,  steady  state 
heat  transfer  analysis  Is  limited  because  It  cannot  account 
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for  solar  gains  and  thermal  mass  capacitance  effects 
Inherent  In  a  structure.  It  lacks  accuracy  In  predicting 
indoor  temperatures  when  other  than  steady  state  conditions 
are  present.  In  order  to  achieve  precise  results  under 
varying  thermal  conditions,  transient  heat  flow  methods  must 
be  applied . 

Transient  Thermal  Conditions 

Transient  thermal  analysis  offers  the  greatest 
precision  possible  for  modeling  heat  transfer  In  buildings. 
The  actual  radiation,  conduction,  and  convection  thermal 
processes  are  accounted  for  with  hour  by  hour  load 
calculations.  Because  of  the  complexity  and  magnitude  of 
the  calculation  process,  computerization  of  this  method  Is  a 
necessl ty . 

Computer  Simulation 

Computer  analysis  Is  a  very  flexible  method  which 
offers  the  user  numerous  options.  Energy  analysis  programs 
can  simulate  heat  transfer  through  an  entire  building 
structure,  and  thus  can  specifically  calculate  the  thermal 
capacitance  effect  due  to  the  building  mass.  Internal  zone 
temperatures  can  be  determined  accurately,  and  hourly 
heating,  ventilating,  and  alr-cond It lonlng  loads  developed. 
Solar  gains  may  be  recorded  from  actual  or  simulated  weather 
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data.  One  of  the  most  prominent  hourly  simulation  computer 
programs  is  called  the  Building  Loads  Analysis  and  System 
Thermodynamics  Program  (BLAST). 

BLAST  Computer  Program 

Developed  by  the  0.  S.  Army  Construction  Engineering 
Research  Laboratory  for  use  by  the  Department  of  Defense, 
BLAST  is  one  of  the  most  widely  used  programs  performing 
transient  thermal  calculations.  It  allows  the  user  to 
describe  the  building  geometry  and  construction  in  detail, 
control  the  input  weather,  and  model  a  number  of  HVAC 
systems.  It  can  simulate  widely  differing  facility  types, 
from  passive  solar  architecture  to  more  typical  office 
occupancy  facilities.  Heat  gain  through  walls,  roof,  and 
slab  is  performed  using  the  ASHRAE  transfer  function  method. 

Although  extremely  capable  for  performing  energy 
analyses,  programs  such  as  BLAST  are  costly  and  time 
consuming.  Computer  simulation  to  design  the  HVAC  equipment 
for  a  multi-zone  office  building  can  be  quite  costly.  This 
expense  is  added  to  the  training  costs  an  engineering  firm 
must  face  before  their  personnel  are  familiar  with  the 
program  input  requirements.  Although  program  execution  is  a 
time  saving  process,  the  time  required  for  training  and  data 
input  prior  to  the  run  can  be  lengthy.  Several  hours  of 
time  would  be  necessary  to  prepare  the  input  deck  for  a 


21 


multi-zone  building.  The  detail  required  of  the  Input 
Information  could  also  be  a  problem.  Programs  of  this  type 
require  hourly  weather  data  which  may  not  be  available  to 
the  user.  The  disadvantages  are  primarily  those  of  cost  and 
time  -  the  benefits  are  the  high  degree  of  accuracy  and 
reliability. 


Method  Comparison 

The  changeover  alternatives  discussed  each  have 
valuable  benefits  for  the  facility  manager  or  energy 
engineer.  Although  the  historical  approach  cannot  determine 
building  thermal  performance,  it  offers  a  system  by  which 
changeover  timing  can  be  evaluated.  When  weather  conditions 
create  thermal  discomfort  for  the  occupants,  the  following 
season's  changeover  date  can  be  adjusted.  Steady  state 
analysis  offers  simplicity  for  the  user  and  a  method  to 
determine  building  temperatures.  Transient  analysis  offers 
precision  and  flexibility.  Solar  and  thermal  mass  gains  are 
calculated  on  an  hourly  basis,  and  the  user  can  input  the 
actual  building  operation  schedule  as  part  of  the  computer 
simulation. 

The  ideal  changeover  method  would  borrow  benefits  from 
each  of  these  approaches.  User  input  information  should  be 
kept  simple  and  to  a  minimum.  Accuracy  should  be  a  high 
priority,  but  it  should  not  depend  upon  expensive  computer 
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analysis  for  calculations.  The  facility  manager  must  be 
able  to  assess  the  thermal  comfort  for  building  occupants, 
yet  be  able  to  maximize  his  savings  by  conserving  energy. 
One  method  to  ensure  thermal  comfort  Is  maintained  Is  to  set 
changeover  timing  as  a  function  of  the  Indoor  building 
temperature . 

Trigger  Temperature  Method 
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be  accomplished  for  the  first  few  hours  of  occupancy,  when 
the  building  Is  the  coldest.  The  resuxtlng  method  combines 
the  accuracy  of  computer  simulation  with  the  simplicity  of 


hand  calculation  procedures.  A  high  precision  level, 
however,  can  only  be  achieved  with  a  carefully  developed 
statistical  analysis  plan. 
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CHAPTER  FOUR 

STATISTICAL  MODEL  DEVELOPMENT 


No  simple  methods  are  available  to  accurately  determine 
Interior  building  temperatures.  The  thermal  capacitance  of 
a  structure  will  cause  a  time  delaying  effect  so  that  gains 
from  solar,  lighting,  and  people  are  not  felt  In  the 
building  Immediately.  Calculations  of  Instantaneous  thermal 
gains  and  losses  can  Inaccurately  represent  the  hourly 
building  temperatures.  Computer  simulation  Is  a  necessary 
tool  to  accurately  describe  the  hourly  temperatures  and 
thermal  conditions.  Through  simulation  of  hypothetical  test 
facilities,  data  may  be  gathered  to  represent  a  broad 
spectrum  of  building  types  and  weather  conditions.  It  Is 
then  possible  to  employ  statistical  analysis  techniques  to 
develop  temperature  predictive  models  which  are  accurate 
within  the  range  for  which  they  are  developed.  These  models 
can  be  formulated  to  retain  precision  without  the  time  and 
expenses  Inherent  In  computer  simulations. 

Statistical  Models 

Three  mathematical  models  may  be  used  to  determine 
building  temperatures:  functional,  control,  and  predictive 
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(Draper  and  Smith,  1981).  If  the  basic  deterministic 
relationships  among  the  variables  can  be  defined,  the 
functional  model  is  best  suited  to  the  problem.  However,  in 
many  cases  where  the  underlying  mathematical  relationships 
are  complex  and  Involved,  the  functional  model  becomes 
difficult  to  develop. 

The  control  model  is  used  in  carefully  designed 
experiments.  As  the  variable  relationships  become  more 
complex,  it  is  difficult  to  "control”  the  variables.  This 
leads  to  the  application  of  the  predictive  model. 

Predictive  models  are  used  when  functional 
relationships  are  complicated  and  difficult  to  control.  For 
example,  they  are  often  used  to  determine  the  strength  and 
reliability  of  construction  materials.  The  compressive 
strength  of  a  concrete  mixture  can  be  calculated  using  a 
functional  model  which  considers  the  chemical  bonding 
processes.  However,  for  the  purposes  required  in 
construction  practice,  the  complexity  Involved  in  this 
procedure  can  be  greatly  simplified  by  using  a  predictive 
model.  An  easier  approach  involves  testing  of  sample 
mixtures  and  recording  the  forces  required  to  produce 
measurable  failures.  From  these  observations,  a  model  is 
developed  to  predict  compressive  strength  based  upon  the 
types  and  quantities  of  the  concrete  mixture  components. 
The  predictive  model  need  not  represent  the  actual  chemical 
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processes  which  cause  concrete  hardening  in  order  to  produce 
an  accurate  and  useful  predictive  model. 

A  predictive  model  can  similarly  be  developed  to 
calculate  building  temperatures  as  a  function  of  the 
envelope,  occupancy,  and  weather  conditions.  The  actual 
thermal  processes.  Including  the  solar  gains  and  capacitance 
effects,  can  be  predicted  as  part  of  the  Indoor  temperature. 
In  order  to  develop  the  predictive  model,  it  is  important  to 
outline  the  analysis  plan  and  describe  the  statistical 
procedures  which  apply. 

Model  Building  Analysis  Plan 

^  An  organized  analysis  plan  is  a  necessary  starting 

i 

point  in  any  problem-solving  procedure.  Methods  to 

establish  goals  and  evaluate  the  final  solution  should  be 
clearly  defined.  The  model  building  process  in  statistical 
studies  can  be  developed  in  three  phases:  planning, 
development,  and  ver if icatlon/ val id  at ion .  Figure  4-1 

presents  the  model  building  flow  diagram  to  be  used  for  this 
statistical  analysis.  The  planning  phase  Includes  defining 
the  problem,  selecting  variables,  outlining  regression 
techniques,  and  establishing  goals. 

Planning  Phase 


The  most  Important  part  in  the  planning  phase  is  the 


Planning 


DevelopBent 


Verification 
and  Validation 
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development  of  a  specific  statement  of  the  problem  (Draper 
and  Smith,  1981).  For  the  work  presented  In  this  thesis, 
the  predictive  model  should  accurately  determine 
transitional  season  building  temperatures  prior  to,  and 
during  the  Initial  hours  of  occupancy.  Changeover  timing 
will  be  based  upon  building  and  weather  conditions  which 
cause  the  building  temperature  to  fall  below  a  user  defined 
threshold.  The  term  "trigger  temperature"  will  be  used  In 
this  thesis  to  describe  the  acceptable  building  temperature 
which  must  be  maintained  for  occupant  comfort.  Above  this 
temperature,  the  building  changeover  can  be  made.  (The 
changeover  condition  has  been  "triggered”.)  Prediction 
models  will  be  developed  on  an  hourly  basis  during  the 
critical  building  hours  from  6-lOam.  The  first  step  In  the 
planning  phase  involves  the  selection  of  possible  model 
variables  . 

Variable  Selection 

The  dependent  variables  for  this  study  are  the  hourly 
building  temperatures.  For  a  facility  functioning  without 
heat,  the  coldest  and  therefore  most  critical  period  Is 
prior  to,  and  during  the  first  few  occupancy  hours.  Four 
dependent  variables  have  been  selected  to  determine  the 
hourly  average  temperatures  from  6-lOam  (Table  4-  1).  This 
enables  the  user  to  assess  building  temperatures  during  the 
critical  occupancy  hours.  | 

i 

I 

I 


29 


Table  4-1 


Model  dependent  variables 

Variable 

Descript  ion 

Units 

TEMP6 

Build Ing 

DBT 

from  6am-7am 

1 

(F) 

TEMP? 

Building 

DBT 

from  7am-8am 

(F) 

TEMPS 

Bui Id Ing 

DBT 

from  8am-9am 

(F) 

TEMP9 

Build Ing 

DBT 

from  9am-10am 

(F) 

Note:  DBT  Indicates  Dry  Bulb  Temperature 

Independent  variables  describe  all  heat  gains  and 
losses  from  occupancy  and  weather  conditions.  These  | 

variables  must  be  fundamental  to  the  problem  and  available  I 

to  the  user  (Draper  and  Smith,  1981).  Building  envelope  | 

variables  should  be  specific  and  In  a  usable  format.  For  j 

example,  Infiltration  can  be  described  In  either  cubic  feet 
per  minute  (CFM)  or  as  air  changes  each  hour.  Both  are 

widely  used  In  load  calculations,  but  the  widely  varying 

range  for  CFM  makes  it  a  poor  choice  for  statistical  I 

i 

prediction  use.  A  large  facility  with  tight  envelope 
construction  could  have  the  same  CFM  Infiltration  rate  as  a 
smaller,  but  older  facility.  The  effect  1,000  CFM  would 
have  on  each  of  these  buildings  would  vary  considerably. 

I 

The  Initial  variable  list,  shown  in  Table  4-2,  Is  developed  j 

from  ASHRAE  standards  widely  used  in  heating  and  cooling 

t 

load  calculations.  ' 


■#p  ■. 


,  T  >-  7* . 
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Table  4-2 

Initial  variables  -  building  and  occupancy 


Variable 

Description 

Units 

APART 

Area  of  partitions 

(sf) 

ASLAB 

Area  of  building  slab 

(  sf  ) 

AWALL 

Opaque  wall  area 

(  sf  ) 

BVOL 

Building  volume 

(Cf) 

OPART 

Density  of  partitions 

( lbs/ c  f ) 

EGLAS 

Projected  area  of  east 

glass 

(sf) 

ELEC 

Lights,  equip  gains 

( watts/ sf ) 

HVAC 

Heat  recovery  system  (70%  eff 

.)  (yes  or  no) 

INFIL 

Infiltration 

(air  changes/hour) 

PSLAB 

Perimeter  of  slab 

(If) 

SCHEDl 

Scheduled  occupancy  begins 

( 7  am  or  Sam) 

SCHE02 

Early  occupancy 

( up 

to  one  hour  early) 

SGLAS 

Projected  area  of  south 

glass 

(sf) 

SHADE 

Glass  shading  coefficient 

( fraction) 

TGLAS 

Area  of  total  glass 

(sf) 

UGLAS 

Glass  U-value 

(Btu/hr*F*sf ) 

UROOF 

Roof  U-value 

(Btu/hr*F*sf ) 

UWALL 

Wa_i  U-value 

( Btu/hr*F*  sf  ) 
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Building  schedule  options  are  best  described  by  a 
special  type  of  Independent  variable,  called  an  Indicator  or 
dummy  variable.  Indicator  variables  are  used  where 
measurements  are  not  recorded  on  a  continuous  and  well- 
defined  scale,  or  when  simplicity  in  model  application  Is 
desired.  The  emphasis  Is  on  qualitative,  rather  than 
quantitative  effects.  Examples  Include  sex  (male  or 
female),  and  political  preference  (Republican,  Democrat,  or 
Independent).  Several  variables  (Table  4-2)  have  been 
described  using  two  level  indicator  variables  (0  *  no 
effect,  1  *  some  effect). 

Weather  conditions  Involve  perhaps  the  most  difficult 
variables  to  select.  A  balance  Is  needed  between  selecting 
sufficient  variables  to  develop  an  accurate  predictive  model 
and  limiting  the  list  to  variables  generally  available  to 
the  user.  Hourly  weather  data  can  provide  very  accurate 
building  energy  calculations,  but  many  potential  users  do 
not  have  access  to  this  Information.  For  this  reason, 
weather  variables  were  selected  on  the  basis  of  daily 
records  (high  and  low  temperatures).  Table  4-3  lists  the 
preliminary  selection  of  weather  variables.  After 
Identifying  the  variables,  statistical  analysis  methods  can 
be  used  to  formulate  the  predictive  models. 
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Table  4-3 

Initial  variables  -  site  and  weather 


Variable 

Description 

Units 

ODBMAX 

Outdoor  dry  bulb 

temp 

max 

(F) 

ODBMIN 

Outdoor  dry  bulb 

temp 

min 

(F) 

SEASON 

Weather  season 

Spring 

or  Autumn 

SOLAR 

Clearness  number 

( f rac  t ion ) 

TGRND 

Ground  temperature 

(F) 

TSDEV 

Temperature  standard 

deviation 

(F) 

Regression  Techniques 


Regression  analysis  Is  the  most  common  statistical 
method  to  predict  responses  from  quantifiable  Independent 
variables  (Draper  and  Smith,  1981).  It  Is  a  widely  used 
method  where  underlying,  but  unknown,  deterministic 
relationships  exist  between  dependent  and  Independent 
variables . 

The  regression  equation,  or  model,  may  contain  one  or 
more  Independent  variables.  The  highest  power  of  any 
independent  variable  is  referred  to  as  the  order  of  the 
model.  The  general  form  of  a  simple,  first-order  model  is: 


Yj  -  Bo  +  BiXj  +  El 


(4.1) 


-  wr  - 


In  this  equation,  is  the  value  at  the  trial,  Bq  and 

are  parameters,  is  the  value  of  the  dependent  variable 
at  the  l^*'  trial,  and  Ej  describes  the  random  error.  The 
error  term  is  the  amount  by  which  any  observation  of  Y  will 
differ  from  the  regression  line.  Linear  regression  models 
are  linear  in  parameters  only  and  may  contain  higher  order 
polynomial  terms.  Because  E  changes  for  each  observation 
and  the  B  terms  can  only  be  determined  after  seeing  all  of 
the  data,  these  terms  are  replaced  with  estimates  in 
regression  analysis.  The  equation  would  then  be  written  as: 

Response  Fit  Error 

Yh  -  b^^Pb^  +  ej  (4.2) 

In  this  equation,  Yj^  is  called  Y-hat,  the  predicted  value  of 
Y,  bQ  and  b^  are  parameter  estimates,  and  X  is  the  value  of 
the  dependent  variable.  The  difference  between  Y  and  Yj^ , 
e^^,  is  the  error  in  the  equation.  This  e^^  term  is  referred 
to  as  the  residual. 

To  obtain  the  best  estimate  of  the  model  parameters, 
regression  analysis  uses  the  method  of  least  squares.  The 
regression  model  which  fits  the  data  best  is  one  in  which 
the  sum  of  squared  error  terms  (SSE)  is  the  lowest.  A 
sample  regression  line  shown  in  Figure  4-2  shows  how  tree 
diameter  may  be  used  to  estimate  the  volume  of  usable  lumber 
which  may  be  produced. 
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In  many  research  problems,  more  than  one  Independent 
variable  Is  required  to  accurately  describe  the  response. 
Multiple  linear  regression  techniques  are  used  for  these 
cases.  The  general  form  of  a  first-order  model  with  two 
Independent  variables  would  be  written: 

Interaction  Term 

Yj  »  Bq  +  +  B2Xj_2  +  ®3*11*12  ®1  (^-3) 

Often  the  response  of  one  Independent  variable  depends 
upon  Its  relationship  with  another  related  variable. 
Interaction  terms  are  used  to  account  for  this  information 
In  multiple  regression  analysis.  In  the  previous  example, 
the  diameter  alone  was  used  to  predict  lumber  volume.  If 
the  Interaction  term,  diameter*helghf  Is  used,  the 
observations  fall  In  a  pattern  closer  to  the  regression  line 
(See  Figures  4-3a  and  4-3b)  .  However,  both  of  these 
patterns  Indicate  nonlinear  trends.  This  can  be  corrected 
by  using  a  variable  transformation.  Transformations 
commonly  used  in  regression  analyses  Include;  1/X^,  1/X, 
l/log(X),  log(X),  and  X  .  The  transformation  Is  clear  In 
this  simple  example.  Figure  4-3c  shows  the  best  fit  of  the 
data  -  a  result  of  using  the  transformation  closest  to  the 
mathematical  equation  for  a  cylindrical  volume.  The 
preliminary  selection  of  interaction  terms  used  In  this 
regression  analysis  is  shown  In  Table  4-4. 
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Table  4-4 

Initial  interaction  variable  terms 


Fred ictor 

Significance 

Interaction 

terms  ■  f(variables) 

Heat  stored  by  partitions 

ADPART 

SB 

APART*DPART 

Infiltration  effect  from 
outside  air  (see  TDIFF) 

AIRCH 

m 

INF1L*TDIFF 

Heat  recovered  from 
early  occupancy  hours 

AIRCHl 

m 

INFIL*SCHED1*HVAC*TDIFF 

Heat  recovered  from 
warm-up  lights 

A1RCH2 

m 

INFIL*SCHED2*HVAC*TDIFF 

Ventilation  effect  from 
heat  recovery  system 

AIRCH3 

m 

INFIL*HVAC*TDIFF 

East  solar  gain  reduced 
by  glass  shading 

EGLASl 

m 

EGLAS*SHADE 

East  solar  effect 

EGLAS2 

m 

EGLAS*SOLAR 

Combined  east  solar 
effect 

EGLAS3 

m 

EGLAS*SHADE*SOLAR 

7am  occupancy  light  gains 

ELECl 

m 

ELEC*SCHEDl 

Warm-up  light  gains 

ELEC2 

m 

::lec*sched2 

Dally  temperature  range 

RANGE 

m 

ODBMAX-ODBMIN 

South  solar  gain  reduced 
by  glass  shading 

SGLASl 

m 

SGLAS*SHADE 

South  solar  effect 

SGLAS2 

m 

SGLAS*SOLAR 

Combined  south  solar 
effect 

SGLAS3 

m 

SGLAS*SHADE*SOLAR 

Losses  through  building 
envelope 

TAVE 

m 

(0DBMIN+0DBMAX)/2 

Balance  point  difference 

TDIFF 

m 

65-ODBMIN 
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Multiple  linear  regression  becomes  a  complex  and 
Involved  calculation  procedure  as  the  number  of  variables 
Increases.  There  are  n*2^  possible  regression  equations  In 
any  analysis,  where  n  Is  the  number  of  Independent  variables 
and  k  Is  the  number  of  dependent  variables  (Neter  and 
Wasserman,  1974).  In  this  study,  the  more  than  thirty 
variable  terms  (variables.  Interactions  and 
transformations),  would  result  In  more  than  a  billion 
possible  regression  equations.  Computer  analysis  Is 
necessary  to  perform  regressions  of  this  scope.  One  of  the 
most  powerful  statistical  packages  available  Is  the 
Statistical  Analysis  System  (SAS).  Developed  by  the  SAS 
Institute  In  North  Carolina,  this  program  Includes  a  variety 
of  statistical  procedures  ranging  from  graphics  to  complex 
multivariate  techniques  (SAS  Institute,  1979).  The 
extensive  descriptive  statistics  of  SAS  make  It  an  Ideal 
tool  for  developing  and  checking  regression  equations 
against  specified  goals. 

Analysis  Goals 

In  order  to  select  the  best  regression  equation  from 
the  many  possible.  It  Is  necessary  to  define  several  goals. 
Preliminary  regression  goals  Include  establishing  the  model 
type  and  maximum  number  of  variable  terms.  A  predictive 
regression  model  will  be  used  in  this  research.  The  model 
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will  be  limited  to  having  twenty  or  fewer  variable  terms 
(parameters)  so  that  hand  calculation  methods  may  be  used  In 
follow-on  applications.  Additional  goals  will  Include 
setting  acceptable  values  for  the  Coefficient  of 
Determination  (R^),  Coefficient  of  Variation,  (C.V.),  and 
parameter  significance  (F  test).  Finally,  models  with 
different  dependent  variables  will  be  compared  for  accuracy 
and  consistency. 

Coefficient  of  Determination 

o 

The  Coeffflclent  of  Determination,  or  R  term.  Is  the 

square  of  the  correlation  coefficient.  It  indicates  what 

proportion  of  the  total  variation  In  the  dependent  variable 

Y  Is  explained  by  the  regression  model  (Walpole  and  Meyers, 

1978).  The  R^  Is  often  given  as  a  fractional  value,  ranging 

2 

from  zero  to  one.  An  R  value  of  zero  means  that  the 
regression  model  Is  of  no  help  In  reducing  the  variation  In 
Y,  and  a  value  of  one  Indicates  a  perfect  fit  with  all 
observations  lying  on  the  regression  line. 

An  R  goal  of  .95  Is  set  to  ensure  that  the  final 
regression  equations  will  be  accurate  predictors  of  the 
building  temperature.  Any  equation  failing  to  meet  this 

minimum  value  can  be  rejected  as  lacking  sufficient 

2 

precision.  One  limitation  of  the  R  statistic  Is  that  as 

2 

the  number  of  pa  ameters  Increases,  the  R  value  cannot 


decrease . 
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Coefficient  of  Variation 

The  Coefficient  of  Variation  is  very  useful  in 
statistical  analyses  because  it  is  a  measure  of  model 
precision  which  takes  into  account  the  number  of  parameters 
in  the  model.  It  is  a  unitless  term  used  to  describe  the 
variation  unaccounted  for  by  a  fitted  regression  equation. 
In  least  squares  regression  analysis,  the  goal  is  to 
minimize  the  sum  of  the  squared  residuals.  One  method  to 
measure  the  variability  of  the  residuals  is  by  calculating 
the  variance  of  the  residuals,  called  the  mean  square  error 
(MSE).  The  Coefficient  of  Variation,  expressed  as  a 
percentage  of  the  mean  response,  is  the  square  root  of  the 
MSE  divided  by  the  predicted  value  mean.  A  value  of  1.0%  is 
set  as  the  goal  for  this  research.  This  will  ensure  that  a 
very  small  amount  of  error  will  occur  in  the  regression 
models.  The  next  step  is  to  ensure  that  all  parameters  are 
statistically  significant. 

Parameter  Significance 

Before  a  regression  equation  can  be  selected  as  the 
best  model.  The  parameters  should  be  tested  for 
significance.  The  possibility  that  a  parameter  is  not  very 
valuable  in  explaining  model  error  is  expressed  In  the  form 
of  a  test  statistic.  First,  a  null  hypothesis  ,  Hq  ,  is 
proposed  that  a  parameter  is  equal  to  zero.  Then  an 
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alternative  hypothesis,  ,  counters  that  the  parameter  Is 
not  equal  to  zero  and  is  indeed  significant  to  the 
regression  equation.  Finally,  an  F  test  statistic  is 
developed  to  see  if  the  parameter  Is  equal  to  zero  within 
specified  probability  limits.  This  test  Is  performed  at 
standard  significance  levels.  The  significance  of  a 
parameter  can  be  evaluated  by  identifying  the  amount  of 
error  which  would  be  added  if  that  parameter  were  dropped 
from  the  model.  In  order  to  ensure  that  the  selected 
regression  models  do  not  contain  any  questionable 
parameters,  the  probability  level  for  the  F  test  is  set  at 
O.OIZ. 

Model  Agreement 

The  final  goal  to  be  met  is  one  of  model  agreement. 
After  regression  equations  for  each  of  the  four  hours  are 
developed,  they  will  be  used  together  to  estimate  the 
building  temperatures  from  6-lOam.  It  is  necessary  that 
these  equations  are  consistent  in  accuracy.  If  one  equation 
differs  from  the  others  substantially  in  its  ability  to 

predict  temperatures,  errors  in  interpretation  could  result. 

2 

The  R  term  is  used  to  measure  accuracy  of  the  final 
predictive  models.  The  variation  in  can  be  used  to 
ensure  that  all  models  are  consistent  with  each  other.  The 
difference  in  the  R  values  for  any  two  predictive  models  is 
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set  at  1.0!S  to  ensure  that  all  models  will  be  consistent  as 
well  as  accurate.  With  all  of  the  goals  set,  the  analysis 
plan  Is  ready  for  the  next  phase  -  development  of  the 
models  . 


Development  Phase 

The  development  phase  Is  central  to  the  statistical 
model  building  process.  During  this  phase,  data  are 
gathered,  and  regression  equations  are  developed,  tested, 
and  evaluated  . 

Data  Collection 

Computer-simulated  data  were  generated  for  94 
hypothetical  single  story  masonry  office  buildings,  using 
the  Building  Loads  Analysis  and  System  Thermodynamics 
(BLAST)  program.  This  program,  developed  by  the  U.  S.  Army 
Construction  Engineering  Research  Laboratory,  uses  the 
ASHRAE  transfer  functions  to  perform  hourly  heat  transfer 
calculations  (U.S.  Army  CERL,  1979).  A  base  condition 
hypothetical  office  facility  of  5000  sf  was  developed  and 
then  modified  through  a  series  of  envelope,  occupancy,  and 
siting  variations.  These  variations  Included  changes  to  the 
building  size,  area  and  U-values  of  envelope  components,  and 
Infiltration  rates,  as  well  as  occupancy  hours  and  site 


orientation. 


A  survey  of  existing  U. 


S. 


Air  Force 


administrative  facilities  was  used  as  a  basis  for  selecting 
the  ranges  for  building  construction  and  occupacy  criteria. 

During  1977  through  1979  the  U.  S.  Air  Force  conducted 
an  Intensl'Te  survey  and  analysis  program  to  reduce  energy 
consumption  In  existing  and  proposed  buildings.  Through 
this  energy  audit  program,  a  data  base  for  Air  Force 
facilities  was  developed.  From  this  data  base  and  Informal 
discussions  with  Air  Force  Civil  Engineering  officials, 
representative  building  areas  and  constructions,  envelope 
component  U-Values,  Infiltration  rates,  and  occupancy 
schedules  have  been  Identified  for  use  In  this  thesis. 
Tables  4-5a  and  4-5b  present  a  listing  of  the  selected 
variable  levels.  A  weather  data  base  was  then  developed  for 
BLAST  input . 

Design  days  for  three  daily  temperature  ranges  (Figure 
4-4)  representative  of  Spring  and  Fall  weather  were 
generated  by  examining  the  29-year  data  base  recorded  by  The 
Pennsylvania  State  University  Department  of  Meteorology. 
The  extreme  high  and  low  recorded  temperatures  were  set  as 
limits  for  the  design  days.  Steps  of  2  F  were  used  to 
develop  58  different  design  days,  having  lows  ranging  from 
18  F  to  54  F  and  highs  ranging  from  28  F  to  84  F. 
Simulation  of  the  hypothetical  test  facilities  was  then 
performed  . 

A  minimum  of  ten  observations  for  each  proposed  model 
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Table  4-5a 

Variable  levels  tested  in  BLAST  simulations 


Levels  -  See  Note  ifl 

Run  #s  Description  (*  Denotes  base  condition) 


1 

to 

8 

Area  of  building 

1800  * 

5000 

12800 

(sf) 

2450 

7200 

16200 

3200 

9800 

20000 

9 

to 

16 

Area  of  partitions 

0  * 

1875 

3750 

(sf) 

625 

2500 

4375 

1250 

3125 

5000 

17 

to 

22 

Density  of  partitions 

6.16 

12 . 5 

32.8 

(Ibs/sf)  * 

12.36 

25.2 

40.7 

59.1 

23 

to 

28 

Area  pec  occupant 

63 

100 

150 

( s  f / per  son) 

75  * 

125 

200 

inf 

29 

to 

35 

Lighting  and  elec  load 

1.5 

2.25 

3.00 

(watts/ sf) 

2.0  * 

2.50 

3.25 

2.75 

3.91 

36 

to 

46 

Site  orientation  angle 

0 

60 

120 

(axis  angle  -  Note  #2) 

30 

*  90 

150 

47 

to 

49 

Occupancy  schedule 

0%  6-7am, 

100% 

7  am-5pm 

(%  lights  and  people)  * 

2  0%  6-7 am, 

100% 

7  am-5pm 

0%  7 -8 am, 

100% 

8  am- 5 pm 

2  0%  7 -8 am. 

100% 

8 am- 5 pm 

50 

to 

53 

Available  sunshine 

.00 

.25 

.75 

(clearness  number) 

.50 

*  .98 

54 

Seasonal  difference 

Autumn 

* 

Spring 

Note  #1  :  The  base  condition  Is  counted  as  a  single 

simulation  and  therefore  Is  not  represented 
as  a  separate  run  with  each  description  group. 

Mote  #2  :  Orientation  levels  were  run  at  20%  and  60% 
glass  levels. 
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Table  4-5b 

Variable  levels  tested  In  BLAST  simulations 

Levels 


Run 

#  s 

Description  ( 

*  Denotes  base 

cond ition) 

55 

to 

61 

Inflltratlon/ventilatlon 

.25 

*  .66 

1.00 

(air  changes/hour) 

.50 

.75 

2.00 

1.00 

3.00 

62 

to 

69 

Heat  recovery  option 

ye  s 

*  no 

(Note  #1) 

70 

to 

74 

Total  glass 

10 

30 

50 

(%  glass  in  ext.  walls) 

*  20 

40 

60 

75 

to 

80 

South-facing  glass 

0 

*  20 

40 

(%  of  south  glass) 

10 

30 

50 

60 

81 

to 

86 

Shading  coefficient 

.24 

.25 

.88 

(Note  #2) 

*  1 .00 

87 

Glass  U-value 

.55 

*1.10 

88 

to 

90 

Wall  U-value 

.0440 

.0786 

*  .1064 

.  2754 

91 

to 

93 

Roof  U-value 

.0409 

.0691 

*  .0897 

.1401 

94 

Base  condition 

Note 

#1 

:  Infiltration  levels  were 

run  at 

30Z  during 

occupied  hours  to  simulate  heat  recovery 
systems . 


Note  #2  :  Shading  coefficients  were  run  for  20X  and 
60Z  glass  levels,  and  for  single  and  double 
glazing  levels. 


Figure  4-4  Daily  temperature  profiles 


variable  is  recommended  for  statistical  regression  analyses 
(Draper  and  Smith,  1981).  For  the  26  Independent  and  4 
dependent  variables,  this  would  require  300  observations. 
Over  1300  observations  were  actually  recorded  so  that  the 
variables  could  be  tested  at  several  levels.  Validation  was 
also  considered  prior  to  data  collection. 

Model  validation  requires  new  observations  to  prove  the 
stability  of  a  developed  regression  equation.  One  of  the 
strongest  and  statistically  useful  validation  procedures  is 
to  use  half  of  the  data  to  generate  the  equation,  and  the 
other  half  for  validation  purposes.  For  this  reason,  over 
2600  observations  were  collected  and  divided  in  half  prior 
to  the  statistical  regression  analysis. 
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I 
I 

I 
1 

Regression  Analysis 

The  regression  model  building  process  Is  a  series  of 
development  and  testing  steps.  After  the  variable  selection 
Is  made,  regression  equations  are  formulated  and  tested 
against  defined  goals.  The  Ideal  regression  equation  is  not 
always  clear,  however.  Often  a  compromise  must  be  made 
between  the  accuracy  desired  and  the  utility  of  the  final 
equation.  The  maximum  number  of  variable  terms  was  limited 
to  20  or  fewer  so  that  the  selected  regression  equations 
would  be  relatively  simple  and  easy  to  use. 

Goals  were  then  set  for  the  statistical  teats.  Any 
models  falling  to  meet  a  minimum  R  value  of  0.95  were 
rejected.  Remaining  models  were  tested  to  meet  a 
Coefficient  of  Variation  of  l.Ot  or  lower.  Parameters  were 
then  tested  for  significance  with  an  F  test  probability  * 

level  of  0.01^.  This  regression  analysis  procedure  was 
followed  for  each  hour  from  6-lOam.  The  resultant  four 
regression  equations  were  finally  compared  for  consistency 
by  making  certain  that  the  R^  values  did  not  vary  by  more 
than  l.OZ  for  any  two  models. 

Appendix  A  Includes  the  Individual  model  statistical 

i 

i 

data  for  all  of  the  generated  predictive  models.  The  | 

detailed  results  of  the  analysis  approach  are  described  In  | 

the  following  sections  for  the  first  hour  dependent 
variable,  TEMP6  (6-7am). 
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Coefficient  of  Detecnlnatlon 

The  first  step  In  testing  the  regression  equation  was 
to  select  only  those  models  with  a  Coefficient  of 

Determination,  or  of  0.95  as  a  mimlnum.  One  of  the  most 

effective  statistical  methods  for  this  purpose  Is  the 
STEPWISE  procedure  available  In  SAS.  Stepwise  regression  is 
a  selective  process  which  chooses  variables  for  model 
Inclusion  based  upon  their  ability  to  produce  a  high  model 
value. 

The  SAS  routine  MAXR  was  selected  to  perform  the 

stepwise  regression  analysis.  This  form  of  stepwise 

regression  begins  with  a  one  variable  model  having  the 

highest  R  value  and  then  selectively  adds  a  second  variable 

2 

to  give  the  greatest  R  Increase.  As  new  variables  are 

added,  MAXR  drops  variables  In  favor  of  more  significant 

ones  so  that  the  "best"  R^  will  be  achieved  for  the  number 

of  variables  present.  Table  4-6  shows  the  results  of 

2 

running  MAXR  for  the  dependent  variable  TEMP6.  The  R 
values  are  listed  by  the  number  of  variable  terms  In  the 
model.  Although  the  stated  goal  of  0.95  was  achelved  with 
only  5  variable  terms,  improvement  Is  seen  until  15  terms 
have  been  Included.  All  equations  meeting  the  R  goal  were 
next  tested  for  error  analysis. 
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Table  4-6  R  Squared  values  for  TEKi-o 


Number  of 
Independent 
Variable 
Terms 

Dropped 

Variable 

Terms 

Ad  d  ed 

Variable  R 

Terms 

Squar 

Value 

1 

— 

TAVE  * 

0,728 

2 

— 

AIRCH  * 

0.883 

3 

— 

SOLAR  * 

0.910 

4 

— 

UGLAS  * 

0.939 

5 

— 

EGLAS3 

0.952 

6 

— 

UWALL 

0.961 

7 

— 

TDIFF 

0.970 

8 

— 

TGLAS 

0.978 

9 

— 

ELECl 

0.983 

10 

— 

UROOF 

0.985 

1 1 

— 

ELEC 

0.987 

12 

— 

AIRCH2 

0.987 

13 

— 

SGLAS2 

0.988 

14 

3GLAS2 

ADPART,  ASLAB^''^ 

0.990 

15 

— 

SGLAS2 

0.991 

16 

— 

INFIL 

0.991 

17 

AIRCH2 

ELEC2,  AIRCHl 

0.991 

18 

ASLABl 

PEOPLE,  ASLAB 

0.991 

19 

— 

ASLAB^^^ 

0.991 

20 

— 

SHADE 

0.991 

Note  * 

Values  failed 

to  meet  goal  of  R^  >  0.95 

Note 
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Coefficient  of  Variation 

The  maximum  value  for  the  Coefficient  of  Variation  was 
set  so  that  the  standard  error  of  the  estimate  would  be  no 
more  than  1.0%  of  the  mean  response.  Two  more  models  were 
eliminated  because  they  failed  to  meet  this  criterion. 
Table  4-7  lists  the  results.  The  remaining  equations  were 
tested  for  parameter  significance. 

Parameter  Significance 

The  parameter  significance  F  test  was  set  not  to  exceed 
a  probability  of  0.01%.  Equations  with  18  and  19  variable 
terms  failed  this  significance  test.  The  remaining 
equations  with  from  7  to  17  variable  terms  met  all 
preliminary  goals  and  were  retained  for  comparison  with 
other  dependent  variable  models.  Summary  statistics  for  all 
models  are  shown  In  Table  4-8. 
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Table  4-7 

Coefficient  of  Variation  for  TEMP6 


Number  of 

Mean 

Standard 

Coeff  of  Var 

(C.V.) 

Ind  ependent 

Square 

Error  of 

C.V.  -  100*3 

/  Mean 

Variable 

Error 

Estimate 

Terms 

(  MSE  ) 

(  s  ) 

(TEMP6  Mean  - 

60.6  F 

5 

0.546 

0.739 

*  1.219 

6 

0.447 

0.669 

*  1 .103 

7 

0. 349 

0.591 

0.975 

8 

0.256 

0.506 

0.835 

9 

0.200 

0.447 

0.738 

10 

0.168 

0.410 

0.676 

11 

0.155 

0.394 

0.650 

12 

0.  144 

0.379 

0.626 

1  3 

0.136 

0.369 

0.609 

14 

0.117 

0.342 

0.564 

15 

0.109 

0.330 

0.545 

16 

0.105 

0. 324 

0.535 

17 

0.102 

0.317 

0.527 

18 

0.  100 

0. 316 

0.522 

19 

0.100 

0. 316 

0.522 

20 

0.099 

0.315 

0.519 

Note  :  * 

Values  failed 

to  meet 

goal  of  C.V.  < 

1  .OX 
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Model  Agreement 

The  last  development  phase  goal  was  one  of  model 
conslsfoncy.  The  goal  was  set  to  have  all  model  R  values 
within  1.02  of  each  other.  Figure  4-5  shows  a  graph  of  the 
R  values  as  a  function  of  the  number  of  variable  terms  In 
the  model.  Plateaus  In  the  TEMP6  regression  equation  can  be 
seen  for  models  with  11  and  15  variable  terms.  In  order  to 
keep  the  TEMP6  R^  value  within  1 . OX  of  the  other  equations, 
the  eleven  variable  term  model  was  selected  as  the  final 
regression  equation. 


Table  4-8 

Summary  statistics  for  prediction  equations 


Number  of  Parameters 


Dependent 

Variable 

Term 

Ind  epend  ent 
Variable 

Te  rms 

R  Squared 
Value 

Coeff  of 
Var  (C.V.) 

Falling 

Significance 

Test 

TEMP6 

11 

0.987 

0.650 

NONE 

TEMP? 

15 

0.979 

0.805 

NONE 

TEMPS 

15 

0.982 

0.694 

NONE 

TEMP9 

16 

0.981 

0.732 

NONE 
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MUICER  OF  VARIABLE  TERHS 


SYMBOLS:  Q)  SELECTED  MODEL 

— -  l.OZ  r2  LIMITS 


Figure  4-5  Model  consistency  comparison 
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The  final  regression  equations  for  all  four  hours  are  then 
written  as  : 

TEMP6  -  (4.4) 

55.25  -.1841*AIRCH  +.2594*TAVE  - . 0851 3*TDIFF 
-.001418*TGLAS+  3.634*SOLAR  -1 5 . 029*UWALL  -20 . 587*UROOF 
-4 . 5821*UGLAS  +.7916*ELEC  +.4407*ELEC1  + . 0085 4 2*EGL AS  3 , 

TEMP7  -  (4.5) 

56.64  -.1734*AIRCH  + . 1 388*AIRCH1  -1.214*INFIL  +.2720*TAVE 
-.1007*TDIFF  +. 01367*EGLAS2  +4.190*SOLAR  -16 . 424*UWALL 
-22 . 431*UROOF  -5.043  *UGLAS  +. 0454*PE0PLE  - . 0000 1 9 2 6* ADPART 
+.5775*ELEC  +2.383*ELECl  - . 000 3305*ASLAB , 

TEMPS  -  (4.6) 

51.07  -.1982*AIRCH  +. 06307*AIRCH1  +. 071 7*AIRCH3  +. 3337*TAVE 
-.02643*RANGE  + . 002 2 37*TGLAS  +. 01 1 0 3*EGLAS2  +4.229*S0LAR 
-16 . 263*UWALL  -25 . 358*UROOF  -3 . 95 37*UGLAS  + . 045 1 9*PE OPLE 
-.00002585*ADPART  +3.102*ELEC  - . 000 3496*ASLAB 
and 

TEMP9  -  (4.7) 

50.51  -.1926*AIRCH  +, 1 353*AIRCH3  +.3353*TAVE 
-.01826*RANGE  +. 004268*TGLAS  +. 01 348*EGLAS2  - . 007 36 7*EGLAS 3 
+5.12*S0LAR  -16 . 769*UWALL  -2 5 . 502*UROOF  -3 . 8904*UGLAS 
+.04750*PE0PLE  - . 0000 30 31 *ADPART  +2.821*ELEC  +.5003*ELEC1 


0004114*ASLAB 
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Verification  and  Validation 

Once  an  equation  has  met  all  of  the  goals  set  in  the 
planning  phase  and  Is  considered  useful,  it  should  be 
examined  for  stability,  analyzed  for  its  practicality,  and 
be  validated  by  an  accepted  method.  The  stability  of  a 
model  depends  upon  the  accuracy  over  its  entire  range  of 
application.  A  model  appearing  to  be  accurate  may  still 
lack  precision  in  a  specific  range.  Verification  that  this 
condition  does  not  exist  is  made  by  examining  plots  of  the 
residuals . 

Realdual  Plots 

Examination  oi  residuals  can  offer  many  Insights  into  a 
selected  model's  stability.  Residual  trends  can  indicate  if 
Important  variables  have  been  omitted,  if  variable 
transformations  are  called  for,  or  If  the  model  variance  is 
not  constant. 

The  residuals,  or  the  differences  between  the 
observations  and  the  predicted  values,  are  plotted  against 
the  dependent  variable.  Figures  4-6a  through  4-6c  show 
residual  patterns  typically  found  In  regression  analyses. 
If  the  residuals  do  not  fall  within  a  horizontal  band,  as  in 
Figure  4-6a,  the  regression  equation  may  require 
modification. 

The  TEMP6  residual  plot  for  the  selected  regression 
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equation  is  shown  In  Figure  4-7.  The  pattern  shows  a 
central  tendency  with  most  points  clustered  between  -1.0  F 
and  +1.0  F.  No  systematic  lack  of  fit  Is  shown  In  the 
residual  pattern.  Residuals  for  the  remaining  dependent 
variables  were  plotted,  with  similar  results  (See  Appendix 
A).  Once  the  prediction  equations  have  satisfied  the 
statistical  requirements,  they  are  analyzed  from  an 
engineering  perspective  to  see  If  underlying  functional 
relationships  have  been  properly  reflected. 

Model  Plausibility 

Engineering  analysis  of  a  prediction  equation  should  be 
made  by  considering  the  model  as  a  whole.  The  least  squares 
regression  coefficients  are  determined  by  adjusting  for  all 
model  variables.  Often  In  multiple  regression  models, 
variables  are  correlated  with  each  other  or  with  variab. 
not  Included  in  the  model.  It  Is  therefore  Important  to 
Interpret  variable  coefficients  from  a  perspective  of  the 
model  as  a  whole.  The  parameters  which  aiw  Included  in  the 
model  should,  however,  be  studied  to  see  If  they  are 
plausible  and  reasonable  predictors. 

Equation  4.4  (page  54)  shows  how  the  building 
temperature  la  a  function  of  envelope,  occupancy,  and  site 
criteria.  The  6-7am  temperature  receives  positive 


temperature  gains  from  solar  radiation  (SOLAR), 


average 
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outdoor  temperature  (TAVE),  east-facing  glass  (EGLAS3)  and 
Internal  electric  loads  (ELEC  and  ELECl).  Losses  are  seen 
In  infiltration  (AIRCH),  the  Indoor  to  outdoor  temperature 
difference  (TDIFF),  the  total  glass  (TGLAS),  and  the 
building  U-Values  (UROOF,  UWALL,  UGLAS).  The  variables 
which  are  not  Included  In  the  selected  model  are  also 
Important  to  consider. 

Although  R  Improvement  Is  seen  in  the  TEMP6  equation 
until  15  variable  terms  have  been  Included  (Table  4-6),  the 
added  terms  do  not  Improve  the  model  significantly.  The 
Infiltration  due  to  a  20%  early  occupancy  option  (AIRCH2) 
and  the  amount  of  south-facing  glass  (SGLAS2)  do  not 
substantially  affect  the  calculation  of  a  building 
temperature.  This  is  possibly  due  to  the  fact  that  the 
early  occupancy  does  not  create  a  significant  infiltration 
load,  and  south  glass  does  not  receive  direct  sun  until 
later  morning  hours.  The  partitions  (ADPART)  are  not 
Important  because  by  6-7am  most  of  the  stored  heat  from  the 
previous  day  has  been  released  to  the  space.  The  building 
area  (ASLAB^^^)  is  also  unimportant  in  improving  the  model 
predictive  ability.  Other  variable  correlations  may  be 
causing  this.  Variable  terms  which  indirectly  describe 
building  area  are  included  in  the  model  before  the  slab  area 
is  considered  (EGLAS3,  TGLAS,  SGLAS2,  ADPART).  They  are 
proportional  to  the  slab  area,  making  ASLAB  '  redundant. 
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All  variables  remaining  In  the  selected  model  are  ones 
developed  from  steady  state  heat  calculation  equations, 
Indicating  that  the  model  la  both  plausible  and  reasonable 
from  an  engineering  standpoint.  Similar  analysis  results 
were  concluded  after  judging  the  engineering  practicality  of 
the  remaining  models.  Before  the  models  can  be  used 
confidently,  however,  they  must  be  validated  using  new  data. 

Validation 

Validation  of  a  statistical  model  Is  a  necessary  step 
In  the  regression  process  (Draper  and  Smith,  1981).  Even 
when  the  selected  equation  is  acceptable  in  statistical 
terms,  Its  ability  to  accurately  predict  future  observations 
should  be  tested.  At  the  outset  of  the  model  building 
process,  half  of  the  data  was  reserved  for  validation 
purposes.  These  observations  are  now  used  to  test  the 
predictive  abilities  of  the  model. 

Two  methods  have  been  used  to  present  the  validation 
summary  for  the  variable  TEMP6.  Both  are  Intended  to 
describe  model  accuracy  In  terms  meaningful  to  the  engineer 
or  end  user  of  the  equation.  The  first  method  defines  the 
model  error  as  a  temperature  deviation  (F)  from  the  BLAST 
results,  and  the  second  method  shows  the  error  as  a  percent 
of  variation.  Figure  4-8  shows  the  frequency  distribution 


for  the  TEMP6  temperature  difference. 


All  observations  for 
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6-7aM  TBMPEiATnRE  DIFFERENCE 


(COMPUTER  -  PREDICTED) 


(  F) 

FREQ 

COM. 

Z 

CUM. 
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z 

-2.0 

1 

0 

0 

0.00 
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1 
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Figure  4-8  TCMF6  temperature  difference  distribution 
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this  hour  are  within  2  F  of  the  simulated  value,  Indicating 
a  high  degree  of  agreement  between  the  predictive  model  and 
the  BLAST  data.  Table  4-9  presents  the  validation  results 
for  all  four  models,  showing  the  distribution  percentages 

Table  4-9 

Summary  of  model  validations 


Description 

TEMP6 

TEMP7 

TEMPS 

TEMP9 

Predictions  within 

1  F  of  BLAST 

97.8% 

91.9% 

94.1% 

9  3.3% 

Predictions  within 

2  F  of  BLAST 

100.0% 

99.4% 

99.5% 

98.8% 

'  Maximum  difference 

(BLAST  -  model) 

-1.9  F 

+2.6  F 

-2.4  F 

-2.6  F 

Average  Error 
(BLAST  -  model)/BLAST 

0.50% 

0.57% 

0.48% 

0.50% 

Maximum  Error 
(BLAST  -  model)/BLAST 

3.1% 

3.7% 

3.4% 

3.6% 

Note  :  Validation  data  contained  1,313  new  observations. 

within  1.0  F  and  2.0  F.  This  table  also  presents  the  second 
method,  showing  the  model  error  as  a  percentage.  The 
average  error  is  less  than  0.60Z  in  all  cases,  and  the 
maximum  error  for  any  observation  is  only  3.7%. 

The  prediction  equations  for  all  models  have  proven  to 
be  precise  and  reliable  predictors  of  building  temperatures 
from  6-lOam.  They  now  can  be  integrated  into  the  Trigger 
Temperature  Method  and  used  in  changeover  analyses. 
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CHAPTER  FIVE 

TRIGGER  TEMPERATURE  METHOD 

The  Trigger  Temperature  Method  Integrates  the  building 
temperature  prediction  equations  developed  in  Chapter  Four 
into  a  system  to  determine  seasonal  changeover  timing.  This 
method  establishes  the  weather  conditions  which  cause 
buildings  to  fall  below  a  minimum  comfort  level,  or  "trigger 
temperature" ,  Average  year  weather  recorded  by  The  Natloal 
Oceanic  and  Atmospheric  Administration  (NOAA)  is  used  to 
predict  the  optimal  changeover  date.  Both  hand  calculation 
and  computerized  methods  are  available  to  the  user.  In  this 
chapter,  two  example  problems  are  discussed  to  demonstrate 
the  Trigger  Temperature  Method  capabilities  and  flexibility. 

Weather  Input  Data 

One  of  the  primary  concerns  in  developing  a  changeover 
analysis  method  is  the  nature  of  the  weather  Information 
required.  In  any  energy  analysis,  weather  data  should  be 
available  and  appropriate  for  the  use  Intended.  Dally 
normal  temperature  data,  developed  by  NOAA  from  30  years  of 
observations,  are  ideal  for  changeover  purposes.  Dally  high 
and  low  temperatures  have  been  recorded  for  numerous 
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locations  and  compiled  Into  representative  "average"  weather 
years.  Table  5-1  shows  the  NOAA  dally  normal  weather  data 
for  Harrisburg,  Pennsylvania,  during  February  through  May. 
(This  table  will  be  used  In  the  example  problems  discussed 
later.)  NOAA  weather  Information  shows  the  seasonal 
progression  on  a  dally  basis,  and  thus  can  be  used  In 
changeover  estimation  methods. 


User  Requirements 


The  Trigger  Temperature  Method  provides  two  procedures 
for  determining  changeover  timing.  If  hand  calculation 
methods  are  required,  linear  programming  techniques  can  be 
used  to  select  the  changeover  date.  If  micro-computer 
capability  exists,  printout  options  are  available  to  provide 
changeover  leclslon  Information.  The  user  flow  chart,  shown 
in  Figure  5-1,  describes  the  Trigger  Temperature  Method  use. 
The  first  step  Is  for  the  user  to  set  temperature  comfort 
goals.  The  Trigger  Temperature  Method  determines  a 
changeover  date  for  a  facility  based  upon  minimum  building 
temperatures  specified  by  the  user.  Trigger  temperatures 
are  set  for  each  hour  to  be  studied.  This  Is  a  minimum  of 
two  hours  -  the  threshold  hour  when  occupancy  begins,  and 
the  upper  limit  when  the  normal  desired  building  temperature 
has  been  reached.  The  threshold  Trigger  Temperature  would 
normally  be  set  at  the  thermostat  night  setback  temperature. 
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Table  5-1 

Dally  normals  of  temperature  for  Harrisburg,  Pa 
(  NOAA  Environmental  Data  Service,  1973  ) 


February 

March 

April 

May 

Day 

Julian 

Julian 

Julian 

Julian 

date 

Max 

Min 

date 

Max 

Min 

date 

Max 

Min 

date 

Max 

Min 

1 

32 

38 

23 

60 

44 

27 

91 

58 

37 

121 

70 

47 

2 

33 

38 

23 

61 

45 

27 

92 

58 

37 

122 

70 

47 

3 

34 

38 

23 

62 

45 

27 

93 

59 

37 

123 

70 

47 

4 

35 

39 

23 

63 

46 

28 

94 

59 

38 

124 

71 

48 

5 

36 

39 

23 

64 

46 

28 

95 

60 

38 

125 

71 

48 

6 

37 

39 

23 

65 

46 

28 

96 

60 

38 

126 

71 

48 

7 

38 

39 

23 

66 

47 

28 

97 

61 

39 

127 

72 

49 

8 

39 

39 

23 

67 

47 

29 

98 

61 

39 

128 

72 

49 

9 

40 

39 

23 

68 

48 

29 

99 

62 

39 

129 

72 

49 

10 

41 

39 

23 

69 

48 

29 

100 

62 

40 

1  30 

73 

50 

11 

42 

39 

23 

70 

48 

30 

101 

62 

40 

131 

73 

50 

12 

43 

40 

23 

71 

49 

30 

102 

63 

40 

132 

73 

50 

1  3 

44 

40 

23 

72 

49 

30 

103 

63 

41 

133 

74 

51 

14 

45 

40 

24 

73 

50 

30 

104 

64 

41 

134 

74 

51 

15 

46 

40 

24 

74 

50 

31 

105 

64 

41 

135 

74 

51 

16 

47 

40 

24 

75 

51 

31 

106 

65 

42 

136 

75 

51 

17 

48 

41 

24 

76 

51 

31 

107 

65 

42 

137 

75 

52 

18 

49 

41 

24 

77 

51 

32 

108 

65 

42 

138 

75 

52 

19 

50 

41 

24 

78 

52 

32 

109 

66 

43 

1  39 

75 

52 

20 

51 

42 

25 

79 

52 

32 

110 

66 

43 

140 

76 

53 

21 

52 

42 

25 

80 

53 

33 

111 

66 

43 

141 

76 

53 

22 

53 

42 

25 

81 

53 

33 

112 

67 

44 

142 

76 

54 

23 

54 

42 

25 

82 

54 

33 

113 

67 

44 

143 

77 

54 

24 

55 

43 

25 

83 

54 

34 

114 

67 

44 

144 

77 

54 

25 

56 

43 

26 

84 

55 

34 

115 

68 

45 

145 

77 

55 

26 

57 

43 

26 

85 

55 

34 

116 

68 

45 

146 

78 

55 

27 

58 

44 

26 

86 

56 

35 

117 

69 

45 

147 

78 

55 

28 

59 

44 

26 

87 

56 

35 

118 

69 

46 

148 

78 

56 

29 

88 

57 

35 

119 

69 

46 

149 

79 

56 

30 

89 

57 

36 

120 

70 

46 

150 

79 

56 

31 

90 

57 

36 

151 

79 

57 
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The  upper  limit  temperature  can  be  set  at  any  value  and  for 
any  hour  after  occupancy  begins.  The  user  also  has  the 
option  of  studying  any  intermediate  occupancy  hour  to  see 
how  the  building  warmup  Is  progressing.  After  determining 
the  hours  to  be  studied,  the  user  selects  the  appropriate 
predictor  equations  and  reviews  the  building  and  site 
characteristics  which  need  to  be  recorded.  The  next  step  is 
the  collection  of  data  required  for  Input  Into  the  predictor 
equations.  Weather  data  for  the  selected  location  are 
gathered  (NOAA  Environmental  Data  Service,  1973),  and  the 
specific  building  characteristics  are  recorded. 

Example  Problems 

Two  example  problems  are  developed  to  demonstrate  the 
different  methods  available  to  the  user.  The  building  Input 
data  for  both  example  problems  are  presented  In  Table  5-2. 

Problem  #1 

The  first  example  uses  linear  programming  techniques 
and  hand  calculations  to  determine  the  Spring  changeover 
date  for  a  hypothetical  5000  sf  office  building  In 
Harrisburg,  Pennsylvania.  Although  this  analysis  Is  made  to 
determine  when  heating  could  be  turned  off  In  the  Spring,  a 
similar  procedure  would  apply  for  turning  the  heat  on  In  the 
Fall.  Occupancy  hours  of  7am-5pm  In  this  example  permit 


1 


68 


Table  5-2 

Building  data  for  example  problems 


Variable  Descrlptlon/Unlt  Value 

APART  Area  of  partitions  (sf)  1875 

ASLAB  Area  of  building  slab  (sf)  5000 

DPART  Density  of  partitions  (Ibs/sf)  12.5 

EGLAS  Projected  area  of  east  glass  (sf)  90 

ELEC  Lights,  equip  gains  (watts/sf)  2.0 

HVAC  Heat  recovery  system  (no)  0 

INFIL  Infiltration  air  changes/hour  .750 

PEOPLE  Number  of  occupants  45 

SCHEDl  Occupancy  hours  (7am-  5pm)  1 

SGLAS  Projected  area  of  south  glass  (sf)  180 

SHADE  Glass  shading  coefficient  0.88 

SOLAR  Clearness  number  (fraction)  .50 

TGLAS  Area  of  total  glass  (sf)  540 

UGLAS  Glass  U-value  ( Btu/hr*F*sf )  0.55 

UROOF  Roof  U-value  ( Btu/hr*F*sf )  .0897 

Wall  U-value  ( Btu/hr*F*sf ) 


UWALL 


.1063 
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analysis  of  four  hours,  6-lOam.  For  this  problem,  trigger 
temperatures  have  been  set  as  follows: 


Trigger  Temperature  Variable 


55 

F 

for 

6-7  am 

TEMP6 

68 

F 

for 

7-8am 

TEMP7 

70 

F 

for 

8-9am 

TEMP8 

72 

F 

for 

9-lOam 

TEMP9 

By  substituting  Into  the  prediction  equations  from  page  54 
and  solving  In  terms  of  ODBMAX  •  f(ODBMIN),  where  ODBMAX  and 
ODBMIN  are  the  outdoor  dry  bulb  high  and  low  daily 
temperatures,  the  prediction  equations  are  reduced  to  ones 
having  with  only  two  variables: 


ODBMAX  •  -2 .721*0DBMIN  +  145.6  (6-7am)  (5.1) 
ODBMAX  -  -2 .696*ODBMIN  +  204.2  (7-8am)  (5.2) 
ODBMAX  -  -2 . 424*0DBMIN  +  192.4  (8-9am)  (5.3) 
ODBMAX  -  -2 . 21 1*0DBMIN  +  186.1  (9-lOam)  (5.4) 


These  equations  can  then  be  graphed,  as  shown  in  Figure 
5-2.  Dashed  lines  show  10  F  and  30  F  daily  temperature 
ranges,  the  limits  for  which  the  prediction  models  were 
validated.  The  Trigger  Temperature  lines  are  plotted  next 
to  determine  the  critical  area  during  which  heat  will  not  be 
required.  (The  critical  area  Is  shown  in  Figure  5-2  as  a 
gray  band.)  The  user  then  selects  possible  changeover  dates 
from  Spring  weather  days  which  fall  In  the  critical  region. 


]  I  I  I  I  •  •  •  I  •  I  I  •  •  •  .  .  .  .  •  I  .  I  I . .  r  ........  .  TTrrrryrrx  .  r.  .TTy 

25  38  38  48  45  58  55  68 

fONIMUM  TEMPERATURE  F  -  SPRIN6 

SYMBOLS  :  ^  HARRISONBURG  DAILY  NORMAL  TEMPS  (NOAA) 

Figure  5-2  Spring  changeover  graph 


Problem  #2 


Although  hand  calculation  procedures  provide  a  quick 
and  simple  method  to  determine  changeover  timing,  the  user 
may  want  to  consider  hourly  building  temperature  profiles 
when  making  changeover  decisions.  This  would  involve 
calculating  the  hourly  building  temperatures  throughout  the 
entire  transitional  season.  Although  this  can  be  done  with 
hand  calculation  or  programmable  calculator  methods,  it  is  a 
tedious  and  time  consuming  process  which  could  be 
conveniently  handled  by  a  micro-computer. 

In  this  example  problem,  the  same  hypothetical  office 
building  has  been  used,  but  the  trigger  temperatures  now 
have  been  set  at  55  F  and  68  F  for  all  hours.  For  the 
computer  application  of  the  Trigger  Temperature  Method,  the 
user  Inputs  the  weather  as  well  as  the  building  data, 
indicating  the  trigger  temperatures  and  applicable  hours  to 
be  studied.  The  summary  output  (Figure  5-J)  selects  the 
weather  days  for  which  the  most  stringent  trigger 
temperature  condition  is  met  and  marks  (*)  the  satisfactory 
building  temperatures.  The  temperature  profile  output 
(Figure  5-A)  graphically  plots  the  hourly  building 
temperatures  during  the  transitional  season  being  studied. 
With  micro-computer  applications,  the  user  is  able  to  select 
a  changeover  date  from  the  hourly  temperature  profiles. 


■UlLDmC  COHPONBIITS  lUILDIMC  AIBAS  BUILDING  CAINS/LOSSES 
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Figure  5-3  Trigger  Temperature  Method  summary  output 


Figure  5-4  Temperature  profile  computer  output 
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Range  of  Application 

As  demonstrated  by  the  two  previous  example  problems, 
the  Trigger  Temperature  Method  provides  a  valuable  tool  for 
facility  managers  or  engineers  to  evaluate  heating  needs 
during  Fall  and  Spring  seasons.  This  methodology  Is  based 
upon  maintaining  a  user  defined  comfort  level  when  the 
building  does  not  benefit  from  an  available  heating  system. 
It  can  provide  accurate  and  useful  results  when  applied 
within  the  conditions  for  which  the  prediction  models  have 
been  validated.  The  user  must  be  careful  to  limit 
applications  of  the  Trigger  Temperature  Method  to  those 
within  the  range  of  applications  listed  In  Table  5-3. 


Table  5-3 


Trigger  Temperature  Method  range  of  application 


INDEPENDENT  VARIABLES 


Oescrlptlon/Unlts 

Minimum 

Maximum 

Building  e 

n  V  e 

lope 

Wall  U-value  (Btu/hr*F*af ) 

0.0440 

0.2754 

Roof  U-value  (Btu/hr*F*sf ) 

0.0409 

0.1401 

Glass  U-value  (Btu/hr*F*8f ) 

0.58 

1.10 

South  glass  (Z  of  south  wall) 

10 

60 

East  glass  (Z  of  east  wall) 

10 

60 

Total  glass  (Z  of  gross  wall  area) 

10 

60 

Glass  shading  coefficient 

0.24 

1.0 

Masonry  construction  (density  Ibs/sf) 

50 

100 

Internal 

Gal 

n  3 

Lights,  equip  gains  (watts/sf) 

1.5 

3.9 

Infiltration  (air  changes/hr) 

0.25 

3.0 

Occupancy  schedule 

7am-5pm 

8am-5  pm 

Occupants  (people/1000  af) 

0 

14 

Area  of  partitions  (sf) 

0 

15700 

Density  of  partitions  (Ibs/sf) 

6.16 

59.  1 

Heat  recovery  system  (  70  Z  eff 

) 

no 

yes 

Weather  F 

act 

0  r  s 

Clearness  number  (fraction) 

0.00 

0.98 

Ground  temperature  (F) 

44 

59 

Outdoor  dry  bulb  temp  max  (F) 

28 

84 

Outdoor  dry  bulb  temp  min  (F) 

18 

54 

Dally  temp  range  (F) 

10 

30 

DEPENDENT  VARIABLES 


Variable  Deacrlptlon/Units  Minimum  Maximum 


TEMP6 

Indoor 

temp 

f  rom 

6-7am 

(F) 

47 

69 

TEMP7 

Indoor 

temp 

f  rom 

7-8am 

(F) 

52 

> 

78 

TEMP8 

Indoor 

temp 

from 

8-9am 

(P) 

54 

> 

78 

TEMP9 

Indoor 

temp 

from 

9-lOam 

(F) 

56 

> 

78 

76 


CHAPTER  SIX 

CONCLUSIONS  AND  RECOMMENDATIONS 


Summary 


Exl 

ating 

me 

thods  whlc 

h  could  be 

used  for  changeover 

estlmatl 

on  h 

ave 

serious 

deficiencies 

in  accuracy 

and 

appl icat 

ion . 

Steady  st 

ate  heating 

load  calculat 

ion 

procedur 

es  1 

ack 

precision 

in  calculat 

ing  hourly  build 

ing 

temperat 

ures 

because  solar 

gains  and  b 

ullding  capaclta 

nee 

effects 

are 

not 

considered 

.  Computer 

simulations  perf 

orm 

accurate 

heat 

transfer  anal 

yses,  but  ar 

e  time  consuming 

and 

expenslv 

e  to 

use  . 

Cha 

ngeov 

er  d 

eclslons  a 

re  usually  t 

he  responsibility 

0  f 

a  faclli 

ty  ma 

nager.  A  me 

thodology  to 

determine  changeo 

ver 

timing  s 

hould 

pr 

ovide  the 

facility  man 

ager  a  procedure 

to 

maximize 
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d  by 

eng Ineer  s 

and  non-engi 

neers  alike. 

ConcluglooB 

The  key  concern  In  a  changeover  method  la  the  level  of 


occupant  thermal 


comfort . 


Changeover  timing 


can  be 
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developed 

In  terms 

o  f 

an  occupant  comfort  threshold. 

or 

" trigger 

temperature”  . 

When  comfort 

conditions  can 

be 

maintained 

without 

mechanical  heating. 

changeover  can 

be 

made  . 

Statistical  analysis  methods  can  be  used  to  determine 
the  weather,  building  envelope,  and  occupant  use 
characteristics  critical  to  facilities  during  transitional 
seasons.  Through  multiple  regression  techniques,  prediction 
equations  can  be  generated  to  calculate  building 
temperatures . 

The  Trigger  Temperature  Method  provides  an  accurate, 
yet  straightforward  method  for  changeover  determination. 
User  Input  Information  Is  tailored  to  use  widely  available 
weather  data  and  engineering  standardized  building  envelope 
descriptions.  Both  hand  calculation  and  micro-computer 
applications  are  offered. 

Recommendations  for  Future  Study 

The  Trigger  Temperature  Method  should  be  extended  to 
Include  different  building  types  and  occupancies.  The  BLAST 
simulated  data  base  should  be  enlarged  to  encompass  varying 
envelope  constructions,  floor  plan  layouts,  and  occupancy 
schedules.  Multiple  regression  analysis  methods  should  be 
used  to  evaluate  new  variables  and  develop  broader,  more 
general  prediction  equations. 


Changeover  criteria  for  air-conditioning  should  be 
developed.  Trigger  temperature  concepts  can  be  extended  to 
describe  cooling  equipment  changeover  timing  in  terms  of 
occupant  comfort. 

Comfort  conditions  for  facility  occupants  should  be 
expanded  to  consider  the  mean  radiant  temperature  effects  in 
a  building.  Trigger  temperatures  can  be  described  in  terms 
of  the  operative  temperatures  to  more  effectively  assess 
occupant  comfort. 

Sensitivity  analyses  should  be  made  to  study  how 
thermal  storage  may  be  designed  to  boost  morning  building 
temperatures.  Thermal  flywheel  effects  could  significantly 
alter  changeover  conditions,  offering  an  extended 
application  for  the  Trigger  Temperature  Method. 

Probability  studies  should  be  performed  on  multi-year 
veather  data.  Trigger  Temperature  prediction  models  should 
be,  used  to  assign  probability  risk  factors  associated  with 
each  calendar  date. 
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APPENDIX  A 

MODEL  DEVELOPMENT  STATISTICS 

The  Trigger  Temperature  Method  uses  statistical  prediction 
models  to  determine  building  temperatures  from  6-lOam.  The 
prediction  models  were  developed  from  multiple  regression 
analysis  of  Building  Loads  Analysis  and  System 
Thermodynamics  (BLAST)  computer  simulations. 

Regression  Statistics 

Regression  models  were  developed  and  tested  against 
goals  defined  for  three  statistics: 

1.  Coefficient  of  Determination  (  R  ) 

2.  Coefficient  of  Variation  (  C.V.  ) 

3.  F  test  for  parameter  significance 

This  appendix  presents  the  statistical  analysis  results 
used  to  select  the  final  prediction  models.  Individual 
model  statistics  and  residual  plots  are  Included.  The 
individual  model  statistics  are  presented  In  Tables  A-1 
through  A-4 .  Residual  plots  are  shown  in  Figures  A-1 
through  A-4.  User  Information  and  variable  descriptions  are 
Included  In  Appendix  B. 
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Table  A-1  TEMP6  regression  statistics 


Number  of 
Independent 
Variable 
Terms 

Dropped 

Variable 

Terms 

Added 

Variable 

Terms 

Variable 

Terms 

Failing 

F  test 

R 

Squared 

Value 

C.  V. 

Value 

1 

TAVE 

0.728 

2.909 

2 

— 

AIRCH 

— 

0.883 

1 . 904 

3 

— 

SOLAR 

— 

0.910 

1.676 

4 

— 

DGLAS 

— 

0.939 

1  .  376 

5 

— -  — 

EGLAS3 

— 

0.952 

1.219 

6 

— 

UWALL 

0.961 

1 .  103 

7 

— 

TDIFF 

— 

0.970 

0.975 

8 

— 

TGLAS 

— 

0.978 

0.835 

9 

— 

ELECl 

— 

0.983 

0.738 

10 

-  — — 

UROOF 

— 

0.985 

0.676 

11 

— 

ELEC 

... 

0.987 

0.650 

12 

— 

AIRCH2 

— 

0.987 

0.626 

1  3 

—  —  — 

SGLAS2 

— 

0.988 

0.609 

14 

14 

SGLAS2 

ADPART 

ASLAB^/2 

— 

0.990 

0.564 

15 

— 

SGLAS2 

— 

0.991 

0.545 

16 

— 

INFIL 

— 

0.991 

0.535 

17 

17 

AIRCH2 

ELEC2 

AIRCHl 

— 

0.991 

0.527 

18 

18 

ASLAB^/^ 

PEOPLE 

ASLAB 

— 

0.991 

0.522 

19 

19 

— 

aslab^/2 

ASLAB 

ASLAB^''^ 

0.991 

0.522 

20 

20 

20 

SHADE 

ASLAB  , 

ASLAB^/^ 

SHADE 

0.991 

0.519 
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Table  A-2  TEMP?  regression  statistics 


Humber  of 
Independent 
Variable 
Terms 

Dropped 

Variable 

Terms 

Added 

Variable 

Terms 

Variable 

Terms 

Failing 

F  test 

R 

Squared 

Value 

C.  V. 

Va  1  ue 

1 

TAVE 

0.487 

3.953 

2 

— 

EGLAS2 

— 

0.673 

3.157 

3 

— 

ELECl 

— 

0.796 

2.493 

4 

— 

AIRCH 

— 

0.889 

1.841 

5 

— 

AIRCHl 

— 

0.921 

1.553 

6 

DGLAS 

0.938 

1 .  382 

7 

— 

SOLAR 

— 

0.956 

1 .166 

8 

— 

UMALL 

— 

0.964 

1.054 

9 

— 

ADPART 

— 

0.969 

0.977 

10 

— 

TDIFF 

— 

0.972 

0.  919 

11 

OROOF 

0.975 

0.875 

12 

— 

PEOPLE 

- - 

0.976 

0.861 

13 

— 

ASLAB 

— 

0.978 

0.826 

14 

— 

INFIL 

— 

0.978 

0.815 

15 

— 

ELEC 

— 

0.979 

0.805 

16 

... 

SGLASl 

0.979 

0.799 

17 

--- 

AIRCH2 

--- 

0.980 

0.795 

18 

SHADE 

18 

SGLASl 

SGLAS2 

0.980 

0. 792 

19 

— 

TGLAS 

TGLAS 

0.980 

0.791 

20 

TGLAS 

20 

EGLAS3 

20 

— 

EGLAS3 

SHADE 

0.980 

0.790 

unclassified 


,  -  INST  OF  TECH  «RIOMT«F*TTlllSON  m  OH  F/0  IS/O 

^JKTAMLOOr  FOO  OCTCRMtHINO  OUILOINO  HEATINO  SYSTEM  CHANOEOVE— ETC(U) 
TTUO  St  ■  a  FORMHALT 

AFIT/el/NR/tl>*OT  NL 
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Table  A-3  TEHP8  regression  statistics 


Number  of 

Variable 

Independent 

Dropped 

Added 

Terms 

R 

Variable 

Variable 

Variable 

Failing 

Squared 

Terms 

Terms 

Terms 

F  test 

Value 

1 

— 

TAVE 

— 

0.517 

2 

— 

EGLAS2 

— 

0.794 

3 

— 

AIRCH 

— 

0.873 

4 

— 

A1RCH3 

— 

0.905 

5 

— 

ELEC 

— 

0.921 

6 

— 

SOLAR 

— 

0.935 

7 

— 

UWALL 

--- 

0.943 

8 

ADPART 

8 

AIRCH3 

AIRCHl 

— 

0.952 

9 

UGLAS 

9 

UWALL 

TGLAS 

— —  — 

0.962 

10 

— 

UWALL 

— 

0.971 

11 

— 

UROOF 

— 

0.975 

12 

AIRCH3 

12 

AIRCHl 

RANGE 

— —  — 

0.977 

13 

— 

AIRCHl 

— 

0.979 

14 

— 

SHADE 

— 

0.980 

15 

PEOPLE 

15 

SHADE 

ASLAB 

— 

0.982 

16 

— 

EGLAS3 

— 

0.983 

17 

— 

INFIL 

— 

0.984 

18 

SHADE 

18 

——— 

SHADE 

EGLAS3 

0.984 

19 

SHADE 

19 

EGLAS3 

19 

ELEC2 

ELEC2 

0.984 

20 

SHADE 

20 

EGLAS3 

20 

ELEC2 

20 

— 

SGLASl 

SGLASl 

0.984 

C.  V. 

Value 

3.571 
2 . 332 
1.835 
1 . 584 
1.447 
1.315 
1 . 225 

1 . 129 

1.002 

0.875 

0.820 

0.785 

0.744 

0.723 

0.694 

0.670 

0.662 

0.660 

0.660 

0.660 
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Table  A'4  TEMP9  regression  statistics 


Number  of 

Variable 

Independent 

Dropped 

Added 

Terms 

R 

Variable 

Variable 

Variable 

Falling 

Squared 

Terms 

Terms 

Terms 

F  test 

Value 

1 

— 

TAVE 

— 

0.468 

2 

— 

TGLAS 

— 

0.782 

3 

— 

AIRCH 

— 

0.841 

4 

— 

SOLAR 

— 

0.873 

5 

— 

AIRCHl 

— 

0.904 

6 

_ 

ADPART 

. — 

0.927 

7 

— 

ELEC 

— 

0.944 

8 

— 

UGLAS 

— 

0.956 

9 

— 

DWALL 

— 

0.965 

10 

— 

EGLAS2 

— 

0.970 

11 

— — — 

DR  OOF 

0.974 

12 

— 

EGLAS3 

— 

0.976 

13 

— 

A1RCH3 

0.977 

14 

RANGE 

14 

AlRCHl 

ELECl 

0.978 

15 

— 

INFIL 

0.979 

16 

ASLAB^^^ 

16 

INFIL 

PEOPLE 

16 

ASLAB^'"^ 

ASLAB 

0.981 

17 

— 

INFIL 

— 

0.982 

18 

— 

ASLAB^''^ 

— 

0.982 

19 

— 

SGLASl 

SGLASl 

0.982 

20 

— — — 

SHADE 

SGLASl 

0.982 

20 

SHADE 

C.  V. 

Value 


3.870 

2.478 

2.116 

1.892 

1.645 

1.439 

1.262 

1.119 

1.003 

0.917 

0.866 

0.828 

0.804 


0.785 

0.779 


0.732 

0.725 

0.719 

0.718 

0.717 


6AH-7AM  PREDICTED  TEMPERATURE  (n 


TEMP6  residual  plot 


r  r  I  T  I  t  I  I  I  I  I  r  r  r  r  I  j  t  rt  TT  rr  r  rfT  fir-t  1  I  I  t  p  r-i  I  I  I  I  I  r  [-r  r  t  i  i  i  i  i  i  j  it  r  rr 


7AM-8AM  PREDICTED  TEHPERATURE  (f) 


TEMP7  residual  plot 


vrr  r  rtt-n 


rtT-f-|-j'rri 


65  70  75 

8AM-9AM  PREDICTED  TEHPERATURE 


TBNP8  residual  plot 


f-r  t-rrr‘|-rT-r-i'T  f-rn  jT'rTTTrT-rrprr-m  t- r  rT'j'n  r  r  rr  rr  t  *  r  rr 
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9AM-I 0AM  PREDICTED  TEMPERATURE  (?) 

TEMP9  residual  plot 
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APPENDIX  B 

USER  INPUT  INFORMATION 

This  apperllx  provides  user  Input  Information  for  the 
Trigger  Temperature  Method.  Descriptions  of  the  variables 
and  the  prediction  equations  for  the  hours  of  6-lOam  are 
Included . 


Variable  Information 

Two  types  of  variables  are  used  In  the  Trigger 
Temperature  Method  prediction  equations.  Continuous 
variables  are  expressed  In  the  engineering  units  shown  In 
Table  B-1 .  Indicator  variables  are  used  to  describe  several 
building  operation  schedule  options.  Tables  B-2  and  B-3 
present  user  Information  for  the  two  indicator  variables 
Included  in  the  prediction  equations. 

Prediction  Equations 

The  statistical  regression  analysis  described  In 
Chapter  Four,  developed  the  final  prediction  equations  as  a 
function  of  variables  and  transformations  (variable  terms) . 
In  order  to  simplify  the  procedures  which  would  be  required 
for  hand  calculation  methods,  all  prediction  equations  have 


been  expressed  In  terms  of  single  variables  only, 
model  prediction  equations  are  listed  In  Table  B-4 . 


Fin 


Table  B-1 

Continuous  Variables 


Variable 

Description 

Units 

APART 

Area  of  partitions 

s  f 

ASLAB 

Area  of  building  slab 

s  f 

OPART 

Density  of  partitions 

Ibs/cf 

EGLAS 

Projected  area  of  east  glass 

sf 

ELEC 

Lights  and  equip  gains 

watts/ sf 

INFIL 

Infiltration  air 

changes/hour 

SGLAS 

Projected  area  of  south  glass 

8  f 

SHADE 

Glass  shading  coefficient 

fraction 

SOLAR 

Clearness  number 

fraction 

TGLAS 

Area  of  total  glass 

s  f 

ODBMAX 

Outdoor  dry  bulb  temp  max 

F 

ODBMIN 

Outdoor  dry  bulb  temp  min 

F 

UGLAS 

Glass  U-value 

Btu/hr*F*sf 

UROOP 

Roof  U-value 

Btu/hr*F*8f 

UWALL 

Wall  U-value 

Btu/hr*F*sf 
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Table  B-2 

Final  model  Indicator  variables 


Variable  Description  (use) 

Value 

HVAC  (Heat  recovery  system  option) 


0  No  heat  recovery  system  available 
1  70Z  efficient  heat  recovery  system 

(Recaptures  heat  by  warming  ventilation  air) 


SCHEOl  (Hours  of  full  occupancy  and  lighting) 


0  8am  to  Spm  -  lOOZ  occupancy  and  lighting  loads 
1  7am  to  Spm  -  lOOZ  occupancy  and  lighting  loads 


Table  B-3 

Building  schedule  options 


Variable 

Schedule 

Solar  time 

06-07 

07-08 

08-17 

V  ax  uc 

:  9 

Opt lo  Q 

SCHEOl 

m 

0 

1 

People 

% 

0 

0 

100 

HVAC 

m 

0 

Lights 

Z 

5 

20 

100 

Recov . 

Z 

0 

0 

0 

SCHEDl 

m 

0 

2 

People 

z 

0 

0 

100 

HVAC 

- 

1 

Lights 

z 

5 

20 

100 

Recov . 

z 

0 

70 

70 

SCHEDl 

m 

1 

3 

People 

z 

0 

100 

100 

HVAC 

■ 

0 

Lights 

z 

20 

100 

100 

Recov  . 

z 

0 

0 

0 

SCHEDl 

m 

1 

4 

People 

z 

0 

100 

100 

HVAC 

m 

1 

Lights 

z 

20 

100 

100 

Recov . 

z 

0 

70 

70 

outo  OUlO  OLnO  OmO 
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Table  B-4 

Trigger  Temperature  prediction  equations 


TEMP6  - 

49.72  +  ODBMIN*( .2148  +.1841*INFIL)  +  .1297*ODBMAX 
-11.97*INFIL  +  S0LAR*(3.634  +.008542*EGLAS*SHADE) 
-.001418*TGLAS  -  20.587*UROOF  -  15.029*UWALL 
-4.5821*UGLAS  +  ELEC*(.7916  +. 4407*SCHEDl ) 


TEMP7  - 

50.10  +  INFIL*('12.48  +9 . 019*HVAC*SCHEDl ) 

+0DBM1N*( .2367  +  INF1L*(.1734  -. 1 388*HVAC*SCHED1) ) 

+. 1360*0DBMAX  +S0LAR*(4.19  +. 01 367*EGLAS )  -22 . 431*UR00F 
-16.424*aWALL  -  5.0432*UGLAS  -  1 . 926*APART*DPART*10**-5 
+ELEC*(.5775  +2 . 383*SCHED1)  -  33. 05*ASLAB*10**-5 
+.0454*PE0PLE 


TEMP8  - 

51.07  +  INFIL*(-12 .88  +HVAC*(4.661  +4 . 010*SCHED1 ) ) 
+0DBMIN*( . 1933  +INFIL*( . 1982  -HVAC*(.07l7  +. 06309*SCHEDl ) ) ) 
+. 1404*0DBMAX  +S0LAR* (4 . 229  + . 01 1 03*EGLAS )  + . 002 2 37*TGLAS 
-25 . 358*UR00F  -16 . 263*UWALL  -3 . 9537*UGLAS  + . 045 1 9*PE0PLE 
-2 . 585*APART*DPART*10**-5  +3. 102*ELEC  -34 . 96*ASLAB* 1 0**-5 


TEMP9  - 

50.51  +  0DBMIN*( . 1859  +INFIL*( . 1926  -. 1 353*HVAC) ) 

+. 1494*0DBMAX  +INFIL*(-12 . 52  +8.794*HVAC)  + . 004268*TGLAS 
+SOLAR*(5.122  +EGLAS* ( . 01 348  - . 007 367*SHADE ) ) 
-25.502*0R00F  -16 . 769*UWALL  -3 . 8904*DGLAS  +. 04750*PEOPLE 
-3.031*APART*DPART*10**-5  +ELEC* (2 . 82 1 0  +. 500 3*SCHED1 ) 
-41 . 14*ASLAB*10**-5 


